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LONG-TERM BEHAVIOR OF REACTION-DIFFUSION EQUATIONS WITH 
NONLOCAL BOUNDARY CONDITIONS ON ROUGH DOMAINS 

CIPRIAN G. GAL AND MAHAMADI WARMA 


Abstract. We investigate the long term behavior in terms of finite dimensional global and expo¬ 
nential attractors, as time goes to infinity, of solutions to a semilinear reaction-diffusion equation on 
non-smooth domains subject to nonlocal Robin boundary conditions, characterized by the presence 
of fractional diffusion on the boundary. Our results are of general character and apply to a large 
class of irregular domains, including domains whose boundary is Holder continuous and domains 
which have fractal-like geometry. In addition to recovering most of the existing results on existence, 
regularity, uniqueness, stability, attractor existence, and dimension, for the well-known reaction- 
diffusion equation in smooth domains, the framework we develop also makes possible a number of 
new results for all diffusion models in other non-smooth settings. 
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1. Introduction 

The mathematical theory for global existence and regularity of solutions to the (scalar) reaction- 
diffusion equation is considered a central problem in understanding models of (non-)degenerate 
reaction-diffusion systems for a variety of applied problems, especially in chemistry and biology. 
It is also essential, for practical applications, to be able to understand, and even predict, the long 
time behavior of the solutions of such systems. It is well-known that the asymptotic behavior 
of solutions to (scalar) reaction-diffusion equations can be well described by invariant attracting 
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sets, and, in particular, by a finite-dimensional global attractor, such that, the dynamics of these 
equations, when restricted to these sets, is effectively described by a finite number of parameters 
(see, e.g., the monographs [ini ESI EH [53]). 

Analytical results for most PDEs in the literature nowadays revolve around the most commonly 
found assumption that the underlying physical space 11 C {N > 2) is smooth enough, and that 
at best, the boundary of H, is of Lipschitz class. But this is barely non-smooth, since a Lipschitz 
boundary has a tangent plane almost everywhere. On the other hand, not much seems to be known 
about partial differential equations (except for some scarce results which we will describe below) 
and their long-time behavior in general, when the physical domain H is actually ’’rough”. This is 
the case, for instance, of domains whose boundary has either a fractal-like geometry or domains 
with cusps which are also frequently used in the applications. Indeed, it cannot be expected that 
objects in the real-world, be they are clouds, trees, snowflakes, blood vessels, etc., will possess 
the structure of smooth manifolds [Hj. One of the main technical difficulties nowadays of dealing 
with ’’bad” domains is the scarcity of Sobolev embedding theorems and interpolation results in 
this general context. In fact, for a general non-smooth domain the usual Sobolev embedding and 
density theorems do not hold [IH] (cf. also Section 2). 

Our main goal in this paper is to develop well-posedness and long-time dynamics results for 
reaction-diffusion equations on domains fl with ’’rough” boundaries, and then subsequently recover 
the existing results of this type for the same models that have been previously obtained in the case 
of domains with smooth boundary cAI. Along these lines, we first establish a number of results for 
scalar reaction-diffusion equations, including results on existence, regularity, uniqueness of weak 
and strong solutions, existence and finite dimensionality of global and exponential attractors, and 
existence of Lyapunov functions. To be more precise, we shall be concerned with diffusion processes 
in ’’rough” domains fl, described by the equation 


dtu — An -|- / (n) = 0 in n X (0, oo), 
subject to the following nonlocal Robin boundary condition 

dyuda + {u + 0^ {u))dpL = 0 on cAI x (0, oo), 

and the initial condition 

u (0) = no in 12. 


( 1 . 1 ) 


( 1 . 2 ) 

(1.3) 

In Eqn. (1.1), / = / (n) plays the role of nonlinear source, not necessarily monotone, and 0^ (n) 
is a certain nonlocal operator characterizing the presence of ’’fractional” diffusion along cffl (see 
Eq. (2.7) below). The normal derivative d^u is understood in the sense of (1.4) specified below, 
a denotes the restriction to 511 of the {N — l)-dimensional Hausdorff measure p is an 

appropriate positive regular Borel measure on dVt. In fact, the regularity assumptions we will 
impose on 511 enter through the measure /r in ( |1.2| ). We will make this more precise in Section]^ 
Since for a ’’rough” domain H the boundary 5H may be so irregular that no normal vector can 
be defined, we will use the following generalized version of a normal derivative in the weak sense 
introduced in |13] . Let /x be again a Borel measure on 511 and let E :!!—)> be a measurable 
function. If there exists a function / G such that 


F-Vifdx= / fipdx+ / ipdpL 


I an 


(1.4) 


for all (f G (5^(11), then we say that p, is the normal measure of F which we denoted by N*{F) := /x. 
If N*{F) exists, then it is unique and dN*{^F) = xjjdN*{F) for all G C^(ll). If n G IL^oc (H) and 
N*(Vu) exists, then we will denote by N(u) := N*(Vu) the generalized normal measure of Vxt. 
The derivative dN{u)/da that we denote by dyU will be called the generalized normal derivative 
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of u. To justify this definition, consider the special case of a bounded domain n C whose 
boundary is Lipschitz continuous, u the outer normal to di} and let a be the (natural) surface 
measure on dQ (in this case, a also coincides with If u G (7^(11) is such that there are 

/ G and g G L^(5n, a) with 


Vu-'Vtpdx= / fipdx+ / QLpda 




idn 


for all If G (7^(11), then dN{u) = gda with g = d^u = dujdv. Throughout the following, without 
any mention, for a bounded arbitrary domain O, we will always mean the identity (1.4) for the 
generalized outer normal derivative of u. 

The interest in analysis and modelling of diffusion processes in bounded domains whose (part 
of the) boundary possess a fractal geometry arises from mathematical physics, and dates back to 
the early 1980’s. The first analytical resnlts aimed at understanding transmission problems, which, 
in electrostatics and magnetostatics, describe heat transfer throngh a fractal-like interface (snch 
as, the snowflake), can be fonnd in |39l I40l ITTl 142] , The type of parabolic problems we consider 
also occur in the field of the so-called “hydraulic fracturing”, a frequently used engineering method 
to increase the flow of oil from a reservoir into a producing oil well (see [TT]; cf. also [36] for a 
related application). Further examples are also provided in the book of Dautray and Lions |22j . 
In all these applications, the mathematical model is usually a linear parabolic boundary valne 
problem involving a transmission condition on a fractal-like interface (layer) which is often a Robin 
boundary condition. The reaction-diffusion equation (1.1) on unbounded fractal domains has also 
I , and in |35| , for bounded fractal domains for which the usual Sobolev-type 


been considered in 

inequalities hold and for which (1.1) is eqnipped with homogeneous Dirichlet boundary conditions 


on dfl. The latter contributions devote their attention mainly to some existence results for some 


special cases of nonlinearities. The motivation to consider (1.1)-(1.3) is also inspired by a wider 


and challenging problems aimed at simulating the diffusion of e.g. medical sprays in the bronchial 
tree |161I17|. In this case, the geometry of the underlying physical domain can be simulated by 
some classes of self-similar ramified domains with a fractal bonndary. Oxygen diffusion between 
the lungs and the circnlatory system takes place only in the last generations of the lung tree, so 


that a reasonable diffusion model may need to involve nonlocal Robin boundary conditions (|1.2) 
the top boundary (the smallest structures), see Section (cf. also O O HI). 


on 


It would be extremely useful if one could give a unified analysis of the reaction-diffusion problem 
0 -( |1.3| ) for a large class of rough domains, including the specified families of ’’fractal” domains 
and/or domains with cusps, using only a minimal number of regularity properties for fl, and then 
use these assumptions about the specific form of the domain, leading to specific models, to derive 
sharp results about existence, regularity and stability of solutions. Then, it is also essential, for 
further practical applications, to show the existence of the global attractor for our general model, 
and then to determine whether the dynamics restricted to this global attractor is finite-dimensional 
or not. Among the first important contributions made to understand the linear problem associated 
with Eqns. (!.!)-(1.2), in general bounded open sets with no essential regnlarity assumptions 

(cf. also 


on cAI, can be found in |59j (cf. also |5Tj for related results). In particular, in |59j it is shown 
that the unique solution of the linear problem is given in terms of a strongly continuous (linear) 
semigroup of contraction operators on (fl), that is order preserving, nonexpansive on L°° (Q), 
and ultracontractive (see Section 2; cf. also |59l Sections 3-5]). The first tool used to derive this 
result is the validity of the following inequality, for arbitrary open sets with finite measure, 


fu 


LN-l (Q) 


< Cn (||Vm||j;^2(q) -I- 11 tt 112,2(9^^0-)) , 


(1.5) 
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which holds for any u G H C'c(n), for some Cq, = C{N,Q,) > 0. This crucial inequality is 

due to Maz’ya [39]. We recall that for Lipschitz domains, the optimal exponent on the left-hand 


side of (1.5) is 2N/{N — 2), while for arbitrary open domains the best optimal exponent is 


2N/ (N — 1), see Section 2.2. The second tool is the notion of relative capacity with respect to 12 
jSli, which is a fundamental tool both in classical analysis and potential theory. Its most common 
property is that it measures small sets more precisely than the usual Lebesgue measure. With 
both these tools at disposal, the local Robin boundary condition (that is, 0^ = 0 in Eqn. (1.2)) 
for the linear heat equation has been investigated in 


I13| (and the references therein) under 
the stronger restriction that 12 possesses the extension property of Sobolev functions. The (linear) 
elliptic system associated with equations (1.1)-(1.2) has also been considered in [8]E]UHlIini[20] , also 
without any essential regularity assumptions on 12. These latter references are mainly concerned 
with the existence of weak solutions for these elliptic systems and several apriori estimates. 

In addition to deriving well-posedness and regularity results for our nonlinear model, perhaps the 
study of the asymptotic behavior is equally as important as it is essential to be able to understand, 
and even predict, the long time behavior of the solutions of our system. One object well-suited to 
describe the long-time behavior is the global attractor. Assuming / G (M) is a function which 
satisfies, for appropriate positive constants. 


Cf |s|^ — Cf < f (s) s < Cf |s|^ + c/, p > 1 


and 


( 1 . 6 ) 


(1.7) 


/ (s) > —Cf, for all s G M, 

resembling the same classical assumptions (see, e.g., naEi]), we prove in Section 3.2 that ( jl.lj )- 
(1.3) generates a nonlinear continuous semigroup {5 (2)}j>o acting on the phase space (12) (see 


Corollary 3.9). This semigroup possesses the global attractor w hich is compact in Lp'{Gl), 

bounded in (12) and has finite fractal dimension (see Corollary 4.8). In addition, we show that 
every unique orbit u{t) = S{t)uo, uq G (12), of the parabolic problem (1.1)-(1.3) ’’instanta¬ 
neously” exhibits an improved regularity both in space and time. More precisely, every such weak 
solution becomes a strong solution (see Section 3.3); hence, Ge,^ contains only strong solutions. 
In fact, the aim of the whole Section 3.1 is to verify the existence of (unique) strong solutions 
which are smooth enough (even with more general assumptions than (1.6)-(1.7)). One of the ad¬ 
vantages of this result is that every weak solution can be approximated by the strong ones and 
the rigorous justification of the regularity results for such solutions becomes immediate. In par¬ 
ticular, exploiting Maz’ya inequality (1.5) we extend the scheme developed initially by Alikakos 
|7] for reaction-diffusion equations, to derive that a particular smoothing property holds for all 


weak solutions of (1.1)-(1.3) with a function / satisfying (1.6)-(1.7). Here we emphasize again that 


the original proofs due to Alikakos [7] rely heavily on the application of the Gagliardo-Nirenberg- 
Sobolev inequality which unfortunately is no longer valid in our context. Besides, in our iteration 
scheme we need to get good control of certain constants in such a way that they do not depend 
on (the bad behavior of) dGl. Furthermore, a key point of this analysis is to employ an appropri¬ 
ate fixed point argument, coupled together with a hidden regularity theorem (Appendix, Theorem 
6.3) for a non-autonomous equation governed by an accretive operator, to deduce smooth strong 
solutions which are differentiable almost everywhere on (0, oo). This regularity is crucial to show, 
for instance, that problem (!.!)-(1.3) has a Lyapunov function on the global attractor Gq^^j. (see 
Lemma 4.3 and Theorem 4.5). Clearly, Gb,/i depends on the choice of boundary conditions 0 
and the measure p on cA2 (see, for instance. Theorems (A) and (B) below). At this point one could 
argue that the long-time behavior of system (1.1 )-(|1.3| is properly described by the global attractor. 


However, it is well-known that the global attractor can present several drawbacks, among which we 
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can mention that it may only attract the trajectories at a slow rate, and that the rate of attraction 
is very difficult, if not impossible, to express in terms of the physical observable quantities (see, 
e.g., [50]). Furthermore, in many situations, the global attractor may not be even observable in 
experiments or in numerical simulations. This can be seen, for instance, for the one-dimensional 
Chaffee-Infante equation 

dtu — dd^u + f {u) = 0 

on the interval (0,1), with cubic nonlinearity / (s) = — s, and non-homogeneous Dirichlet bound¬ 

ary conditions (i.e., u{0,t) = u{l,t) = —1, t > 0), in which case every trajectory is exponentially 
attracted to the ’’single point” attractor {—1}. On the other hand, this problem possesses many 
interesting metastable “almost stationary” equilibria which live up to a time ~ e'’* ^ and, thus, 
for d > 0 small, one will not see the global attractor in numerical simulations. This is known to 
happen, for instance, for some models of one-dimensional Burgers equations and models of pattern 
formation in chemotaxis (see |50|, for further references). Henceforth, in some situations, the global 
attractor may fail to capture important transient behaviors. Besides, in the general setting of ar¬ 
bitrary open domains, this feature can be further amplified for our system due to the boundary 
condition (1.2) and the ’’rough” nature of the domain near dQ. It is thus also important to 


construct larger objects which contain the global attractor, attract the trajectories at a fast (typi¬ 
cally, exponential) rate which can be computed explicitly in terms of the physical parameters, and 
are still finite dimensional. A natural object is the exponential attractor (see, e.g., (50] Section 3]; 
cf. also below). In Section we prove the existence of such an exponential attractor not only for 
the dynamical systems generated by the weak solutions (see Theorem 4.7) but also by the strong 
solutions, which require essentially no growth assumptions on the nonlinearity as |s| —)• oo (see 


Theorem 4.II). Roughly speaking, the assumption on / reads: 


lim inf 

|s|—^- 1-00 


f{s) 


> -A* 


( 1 . 8 ) 


for some A* G [0,(70), where Cq = C {N,Q) > 0 is the best Sobolev/Poincare constant into (1.5). 
The latter result seems to be new for (!.!)-(1.3) even when is a smooth domain. We refer the 


reader to Section for the precise assumptions, related results and generalizations. 

We emphasize that the measure /r and/or boundary regularity assumptions we employ are of 
general character, and as a result do not require any specific form of the domain H; this abstraction 


allows (I.I)-(I.3) to recover all of the existing diffusion models that have been previously studied 


in smooth bounded domains (including Lipschitz domains), as well as to represent a much larger 
family of mode ls for (I.I)- ([Tl that have not been explicitly studied in detail. As a result, the 
system in (!.!)-(1.3) includes reaction-diffusion models, on domains Q which possess the 


extension property of Sobolev functions, and non-Lipschitz domains whose boundary is only Holder 
continuous, as special cases, and on many arbitrary open bounded domains H (satisfying, for 
instance, (1.5)) not previously identified. In Section]^ we discuss how the unified analysis presented 
here can be used to establish the same results for other important classes of partial differential 
equations, such as, reaction-diffusion systems for a vector-valued function it = (tti, ...,Uk) ■ 

The following theorems can be treated as special cases of our results. They apply for instance 
to domains H with a ’’fractal” boundary dVl. 

The o rem (A). Assume that H has the -extension property of Sobolev functions and f satisfies 
(1.6)-{1.7). Let IX he the restriction to dLl of the d-dimensional Hausdorff measure for any 
d G {N — 2, N). Then, the dynamical system S {t) : (fl) —)• uq u{t) = S{t)uo, 

associated with weak solutions for the parabolic problem {73-(73) is gradient-like and possesses 
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the global attractor of finite fractal dimension. Moreover, the semigroup S {t) also has an 
exponential attractor 

The orem (B). Assume that Q has the -extension property of Sobolev functions and f obeys 
(l.S). Let pL be the restriction to dLl of the d-dimensional Hausdorff measure for any d G 
{N — 2,N). Then, the dynamical system T (t) : (Q) —)■ L°° (O), uq ^ u{t) = T (t) uq, associ¬ 

ated with ’’strong” solutions of the parabolic problem (1.1\ )-{1.3) has an exponential attractor 
(hence, also a global attractor of finite fractal dimension). 

The assumption on the measure n deserves some additional comments. First, we note that when 
d = N — 1 but jjL = TL^q^ (= a) is locally infinite, that is, for all x G dO. and 5 > 0 we have that 

a [Bs {x) r\ dLl) = oo, ( 1 - 9 ) 


in this case the boundary value problem (1.1)-(1.2) coincides with the (homogeneous) Dirichlet 
problem for (1.1). If (1.9) only holds on part of the boundary Too C dLl, Dirichlet boundary 
conditions are satisfied on that part Fqo and the usual (nonlocal) Robin boundary condition (1.2) 


is satisfied on the remaining part (9D\roo. On the other hand, for domains with fractal boundaries 
dLl, one often finds some do G {N — 2, N) such that 

Ldo = Ldo ( 9 ^) < oo- 

In this case, one may think of ( |1.2[) as a generalized (nonlocal) Robin boundary condition {p, = 
for the boundary value problem (jl.ip-( 1.3). We once again refer the reader to Section]^ for more 


precise statements of the results on the existence of finite dimensional global and exponential 
attractors. 

The remainder of the paper is structured as follows. In Section we establish our notations 
and give some basic preliminary results for the operators and spaces appearing in the model (1.1)- 
(1.3). In Section]^ we prove some well-posedness results for this model; in particular, we establish 
existence results for strong solutions (Section 3.1), weak solutions (Section 3.2) and then derive 
regularity (Section 3.3) and stability results (Section 3.2). In Section Jdjwe prove results which 
establish the existence of global and exponential attractors for 0- ( |1.3| ), and the existence of 
Lyapunov functions. Section contains some concluding remarks. To make the paper reasonably 
self-contained, in Appendix, we develop some supporting material on regularity results for abstract 
non-homogeneous evolution equations which are necessary to derive the results in Section 3.1. 


2. Preliminaries 

2.1. Some facts from measure theory. As mentioned in the introduction, we aim to analyze 
the long-time behavior of solutions to the reaction-diffusion equation (1.1) (supplemented with the 
boundary and initial conditions (1.2)-(1.3)) in terms of global and exponential attractors [50]. In 
order to make the paper as self-contained as possible, in this section we recall the main definitions 
and results in measure and Sobolev function theory which will be extensively used in what follows. 


Let Ll C be an open set with boundary dLl. We say that a Borel measure p on dLl is regular 
if for every A C dLl there exists a Borel set B such that A C B and p{A) = p{B). A measure 
p on DLL is a Radon measure if p is Borel regular and p{K) < oo for each compact set K C dLl. 
Let a be the restriction to dLl of the {N — l)-dimensional Hausdorff measure Then it is 

well-known that cr is a regular Borel measure on dLl but it is not always a Radon measure, that is, 
if D has a “bad” boundary dLl, it may happen that compact subsets of dLl have infinite u-measure 
(see [25l[26]). 
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Definition 2.1. For a regular Borel measure /i on we denote by 

•= {z ^ do. ■. r) n cffi) = cx) V r > 0} (2.1) 

the relatively closed subset of dQ on which the measure g, is locally infinite. The complement of 
r“ denoted by 

:= dVl \ r“ = {z e 3 r > O : p,{B{z, r) n dFt) < 00 } (2.2) 

is the relatively open subset of dQ on which the measure pL is locally finite. We will also call the 
part the domain of the measure p and denote D{p) = 

We notice that the restriction of p to is a Radon measure and p{T^) < 00 if is such that 
its boundary is a compact set. For more properties of Borel measures and the Hausdorff measure, 
we refer the reader to the monographs [251 [26] . 

Example 2.2. It is clear that by definition, if p{d^) < 00 , then D{p) = so that F^ = 0 . 

Next, let 0 < On+i < bn+i < 0 ^ < 1 (n G N) be such that lim„_,.oo an = 0. Let 

^ := I (a;, y) G (0,1) X (0,1) \ y [an, bn] x 
I nSN 

Let a := Td^ldn. It is easy to see that a{dQ,)= 00 and that F^ = {0} x [1/2,1]. 

Example 2.3. Let C be the Koch snowflake (see Figure 1) and let a := 'H^\dn- Then 

F“ = dQ so that F^- = D{a) = 0 . It is well-known that dit has Hausdorff dimension equal 
to d := ln(4)/ln(3), see [55|. In this case, the d-dimensional Hausdorff measure of dQ is finite, 
n'^idn) < 00 , so that = 0 and D{n'^) = F^d = dQ. 

Next, we give another example which is also meaningful for the applications. 

Example 2.4. Let now 11 C be the class of tree-shaped domains with self-similar fractal 
boundary introduced in |l5| . Such a geometry can be seen as a two-dimensional idealization of the 
bronchial tree (see [I] ; cf. also laiaiiiEi). In order to describe this set, we follow [T|. Consider four 
real numbers a, a, jd, 9 such that 0 < a < a > 0, /3 > 0 and 0 < 0 < 7r/2. Let /*, i = 1, 2, be 

two similitudes in given by 

/i (xi, X 2 ) = (—a, fd) + a {xi cos 6 — X 2 sin0, xi sin 6 + X 2 cos 9 ), 

/2 (xi, X 2 ) = (a, (d) + a (xi cos 9 + X 2 sin0, —xi sin 0 -|- X 2 cos 9) . 

One can define by Fj the self-similar set associated to the similitudes fi, f 2 , i-e., the unique compact 
subset of such that Tj = fi {Tj) L) f 2 (Fj). To construct the ramified domain whose boundary 


1 

2 ’^ 


C 
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Figure 2. Tree-domain (0 < a < a*) 

is Fj, we further let An be the set containing all the 2” mappings from n} to {1,2}, called 

strings of length n for n > 1, and define A = |J A„ as the set containing the empty string and 

n>l 

all the finite strings. Next, consider two points Pi = (—1,0) and P 2 = (1,0) and let F;, be the line 
segment [Pi,P 2 ]. We impose that /2 (Pi) and /2 (P 2 ) have positive coordinates (i.e., a cos 9 < a 
and asin0 < j3). The first cell Yq of the tree domain (i.e., the bottom of the tree, see Figure 2) is 
constructed by assuming that Yq is the convex, hexagonal, open domain inside the closed polygonal 
line joining the points Pi, P 2 , /2 (P 2 ), /2 (-Pi), /i (P 2 ), /i (Pi), and Pi in this order. With the 
above assumptions on 0,a, /3, a, this is true if and only if (a — 1) sin 6 + (3 cos 0 > 0. Under these 
assumptions, the domain Yq is contained in the half-plane X 2 > 0 and is symmetric with respect to 
the vertical axis xi = 0. Next, we introduce Kq = Yq. It is possible to glue together Kq, fi {Kq) 
and /2 {Kq) and obtain a new polygonal domain, also symmetric with respect to the axis xi = 0. 
We define the ramified open domain 12, 

n = Interior U fAKo)] , 

VfreA / 

see Figure 2. Note that 12 is symmetric with respect to the axis xi = 0. In [35] it was proved 
that, for any 0 < 9 < tt/2, there exists a unique positive number a* < l/\/2, which does not 
depend on (a,/3) , such that for 0 < a < a*, Fj has ”no self-contact”. In this case the Hausdorff 
dimension of Fj is d = — In (2) / In (a) > 1 and by |5|, possesses the W^’^-extension property of 
Sobolev functions. Moreover, since P'^(Fj) < 00 we have F“^ = 0 and = F-^d = dO, D I'f 

(cf. Definition |2.l[ ). On the other hand, note that for a = we have F“ = Tj so that 

F^ = D{a) = dtr\Tf. 

Definition 2.5. Let d G [0,N] and let fi be a regular Borel measure on dH. 

(a) We say that fi is an upper d-Ahlfors measure if there exists a constant C 2 > 0 such that for 
every x G dLl and every r G (0,1], one has 

IJ,{B{x, r) n dLl) < C 2 r'^. 

(b) We say that fi is a lower d-Ahlfors measure if there exists a constant Ci > 0 such that for 
every x G dO. and every r G (0,1], one has 

/i(P(x, r) n dLl) > Cir'^. 

By [HI Remark 2.20], every upper d-Ahlfors measure is absolutely continuous with respect to 
the d-dimensional Hausdorff measure that is, there exists a function / G L°° such that 

/i(A) = fj^f (x) dlA for every Borel set A C dVt. Moreover, by |43[ Lemma 1.17], the d-dimensional 
Hausdorff measure is absolutely continuous with respect to every lower d-Ahlfors measure /i. 
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Definition 2.6. Let F be a closed set in We say that F is a d-set for some d G (0, A^] in the 
sense that there exist a Borel measure ^ on F and some eonstants C 2 > Ci > 0 such that 

< //(F n B{x, r)) < C 2 r^ for all x € F and all r G (0,1], (2-3) 


where B{x,r) is the Euclidean metric ball. 
Ahlfors measure. 


Notice that a measure satisfying (2.3) is called a d- 


Let now Ll C be an arbitrary bounded domain with boundary dLl. By [38] (see also m), 
dLl is a d-set if and only if (2.3) holds with /r being the restriction to dLl of the d-dimensional 
Hausdorff measure T-L^. Hence, it follows directly from (2.3) that is an upper d-Ahlfors measure. 


In particular, D{l-r) = = dLl so that 

B“ = cAI, where we recall that a = 


= 0. If d > A" — 1, then D{(t) = B^ = 0 so that 


be the K och s nowflake and recall that d = ln(4)/ln(3) is the Haus- 
Example 2.3). By |55|, the restriction of to dVL is an upper 


C be the tree-domain in Example 2.4 By |3Z|, there ex- 


Example 2.7. Let Ll C 
dorff dimension of dLl (cf. 
d-Ahlfors measure. Let now 
ists a unique regular Borel invariant measure /U on Bj in the sense that for every Borel set B, 
/i (B) = n [fi^ (B)) -|- n (/ 2 "^ (B)) . The set B f, endowed with the self-similar measure //, is a d-set 
where d = — In (2) /ln(o) (see [T]), hence /i|p^ = 

Eor more properties of Ahlfors measures and d-sets we refer the reader to Ha EU ESI EH] and 
the references therein. 


2.2. The Maz’ya space. In this subsection, we introduce the function spaces needed to study our 
problem. We recall that Ll C {N > 2) is an arbitrary bounded domain with boundary dQ and 
the measure a = 


Eor p G [1, 00 ), we denote by IE^’P(II) the first order Sobolev space endowed with the norm 

l|w||ivi.p(D) = . 

Then W^’^{Ll) is a Banach space and VB^’^(n) is a Hilbert space. We let 

W1'2(B1) := IEL2(Bl)nC'(I2) 

Then is a proper closed subspace of IB^’^(n) such that they coincide if, for instance, 11 has 

the segment property (see jlH Theorem 1, p.28]). We recall that the latter property is equivalent 
to n being of class C, which is slightly weaker than the Lipschitz property of 11. 

The following important inequality is due to Maz’ya [Ml Section 3.6, p.l89]. 


Theorem 2.8 (Maz’ya). There is a constant C = C{N) > 0 such that for every function u G 

VBi’i(ii)nC(n), 

||w||^^n ^ (l|V“lli,n + ll«llLi(an,<T)) • (2-4) 

In the inequality (2.4), the constant C{N) is exactly the so called isoperimetric constant. 


The following result is a direct consequence of the Maz’ya inequality (2.4). 

Corollary 2.9. f |49l Corollary 2.11.2]j Let 1 < p < 00 and p* :=pN/{N — l). There is a constant 
C = C{N,p, n) > 0 such that for every u G IT^’^(O) n (7(11), 


\u\ 




< C 


p,n 


+ "w 


LP{dn,(7) 


(2.5) 
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Remark 2.10. The exponent p* in (2.5) is optimal in the sense that it cannot be improved without 
further regularity assumptions on 12, see [491 Example 2.11, p.l23]. 


Now, we let the classical Maz’ya space ITp ^(12, dQ) to be the abstract completion of 


W^,p := {u G W^’P{n) n C(f2) : / \u\Pda < oo} 

Jdn 


with respect to the norm 


(o,ao) 


P -L Ib/ll^ 
p,n ^ ll“llLP(an,cr) 


i/p 


By (2.5), it follows that 

VEpp(12, 912) LP (12) (only continuous embedding). 


( 2 . 6 ) 


In this paper, we focus on the Hilbert space case p = 2, that is, the space IT 2 ^ 2 (^) 912). 

Next, let fjL be an arbitrary regular Borel measure on the boundary 912. We define the Maz’ya 
type space W 22 {^-:d^^ d) to be the abstract completion of 


iy^ 2 := {u G ITb2(12)nC(12) ; / |u|^d/x < oo} 


/an 


with respect to the norm 


\wl^{n,an,p) 


= u 


WL2(n) 


+ 


/an 


|u|^ dp 


1/2 


It follows from (2.5) that the spaces 1^2 2(12, 912, cj) and IE 2 ^ 2 (^) ^112) coincide with equivalent norms. 
Now, for a measurable function u on 912 and s G (0,1), we let 

Ksix,y)\u{x)-u{y)fdp^dp , Ks{x,y) := 


_7/1-^—lH-2s 


\x - y\ 


J Jauxdn 

We let the fractional order Sobolev-Besov type space 

H^{dQ,p) := {u G L^(912,/i) : Ms,p{u) < ooj 
be equipped with the norm 

1 /2 

ll'“llL/‘>(an,/.j) := (l|//|lL2(an,At) +-^s,piu)^ 

For further properties of the space R*(912, p) we refer the reader to |21j and the references therein. 
Let (H^(912,//))* denote the dual of the Hilbert space H^{dQ,p). We dehne on the boundary a 
(nonlocal) operator 0^ : H^{dQ,p) —> p))* as follows: for u, a G R®(912,/i) we set 

(0^(u),a) := / / Ks{x,y){u{x) - u{y)){v{x) - v{y))dp^dp (2.7) 

J Janxdn 

where (•, •) denotes the duality between H^{dQ,p) and (77^(912,/i))*. With some further abuse of 
notation, from now on (•, •) denotes the duality between any Banach space X and its dual X*. 
Now, we let V 0 (12, 912,/i) be the abstract completion of 


Wp,e ■■= {u G Wp,2 : 


Ksix,y)\u{x) -u{y)\^dp^dp < 00 } 


anxdn 
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with respect to the norm 


; 1 . 2 / 


+ / + 
Jan 


anxdn 


Ksix,y)\u{x) - 


1/2 


It is clear that V 0 is continuously embedded into VF 2 ^ 2 (^) A*) if /t = cr, then 

\\\u\\\ := (\\S/u\\l^+ [ \u\‘^da+[ [ Ks{x,y)\u{x) - u{y)\‘^daxday') 

V Jan J Jauxan J 

defines an equivalent norm on the space V 0 ^(n, 50, cr), and this space is continuously embedded 
into W 22 {^i^^)- We notice that if fl has Lipschitz boundary, then the spaces = 

11^2 ~ with equivalent norms. Moreover, if the measure /i = 0, then we always 

have VF 2 ^ 2 (^) 0) = V 0 ^(II, 90,0) = >V^’^(n) with equivalent norms. On the other hand, for 

arbitrary domains or/and for arbitrary regular Borel measures, in order to have an explicit descrip¬ 
tion of the Maz’ya type space VF 2 ^ 2(^5 d) and its subspace Vq (0, dVt, y), we need the following 
notion of capacity. 

Definition 2.11. The relative capacity Cap 2 ,Q. with respect to fl is defined for sets A G Q by 

u G VV’^’^(n), 3 0c open, 

3^2) ■ ^ c O and u > 1 a.e. on id GO 


Cap 2 ,n{A) := inf < ||m||^i ,2 


• A set P Gil is called 00112 ,D-polar if Cap 2 ,n{P) = 0. 

• We say that a property holds Cap 2 ,n-Quasi everywhere (briefly q.e.) on a set A Gil, if there 
exists a Cap 2 ,n-polar set P such that the property holds for all x G A \ P. 

• A function u is called Cap 2 ,Q-quasi continuous on a set A C if for all e > 0 there exists 

an open set O in the metric space such that Cap2,n{0) < e and u restricted to A\0 is 
continuous. 

The relative capacity Cap 2 ,n has been introduced in [8] (see also [aiMllSH] for a more general 
version) to study the Laplace operator with linear (local) Robin boundary conditions on arbitrary 
open subsets in . Note that if then Cap 2 ^N = Cap 2 is the classical Wiener capacity 

HUES]. By jSl Section 2], for every u G W^’^(n) there exists a unique (up to a Oop 2 ,n-polar set) 
Oap 2 ,r 2 -quasi continuous function ft : fl —)• M such that u = u a.e. on il. Moreover, if tt G W^’^(n) 

and Un G >V^’^(n) is a sequence which converges to u in then there is a subsequence of 

Un which converges to u q.e. on il. 

We have the following situation regarding the Maz’ya type space 11^2^ 2 h)- Let r“ be the 
relatively closed, and let D{y) = T^ be the relatively open subsets of dil defined in (2.1) and (2.2), 

respectively. Then every function u G is such that ri|r“ = 0, where we recall that 


Wyy := {u G W^’\il)nC{il) : 


\u\ 


Ian 


^"dy < 00 }. 


Let T C dil be a relatively closed set. Since the closure of the set {u G W^’'^{il) H (7(11) : u|r = 0} 
in W^’‘^{il) is the space {u G W^’^(n) : u = 0 q.e. on T}, it follows that every function u in 
IT 2 ^ 2 (^) A*) (see [H] |9]). With these observations, one has the following 

descriptions: 

W2^,2(^!A^) = ^ Ly^’^(O), ft = 0 q.e. on r)f and f \u\^dy<oo} 


(2.8) 
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and 


Ks{x, y)\u{x) - u{y)\^dy.^dix < cx)}. (2.9) 


Note that if r“ = dO. or Cap 2 ^n(r^) = 0, then VF2^2(r2, 90,/x) = V0^(O,9O,/x) = VFq’^(O) := 
-VP1'2(£7) 

'D(O) . This is the case for the well-known 2-dimensional open set bounded by the Koch 

snowflake if one takes the measure y = a and is also the case for d-sets if /x = a and d > Ai — 1. 


Throughout the remainder, we identify each function u G W^’^(O) with its quasi-continuous 
version u. With this identification, we define the bilinear symmetric form with domain 

V0’^(O,9O,/x) by 


Ae^^{u,v) := / Vu-Vndx-I- / uvdy 

Jn Jdn 

Ks{x, y){u{x) - u{y)){v{x) - v{y))di^^dfiy. 


+ 


anxdn 


( 2 . 10 ) 


Since every function u G Vq^(0, 90, y) is such that u = 0 q.e. on T)^, one has that the form 
defined in (2.10) is given by 

Aq^^{u,v) ■.= / Vu-Vndx-I- / uvdy 

J n 'JT u 


+ 

'2 


Ks{x, y){u{x) - u{y)){v{x) - v{y))dn^dy, . 


( 2 . 11 ) 


The form is closed on L"^(0) if and only if the functional ■ T^(0) —)• [0, -|-oo] defined by 




_ J ^Ae,,j,{u,u) if u G Ve’^(0,90,^), 


- 1-00 


if u G T^(0) \ Vq^( 0, 90,/x) 


is lower semi-continuous on L^(0). To characterize the lower semi-continuity of and hence, the 
closedness of ^©,^, we need to following notion of boundary regularity and/or measure regularity. 


Definition 2.12. We say that a subset T o/90 is Cap 2 ,n'^dmissible with respect to /x, if for every 
Borel set ^ C T, one has Cap 2 n(^) = 0 implies n{A) = 0. 

The following result is taken from |59l Theorems 4.4 and 5.2] (see also |56l Theorem 4.2.1] and 
m Theorem 2.11]). 


Theorem 2.13. The following assertions are equivalent. 

(i) The operator R : W 2 2 {^,d^l, y) —)• Lf‘{TL), n I—)> xx|q is injective. 

(ii) The set T^ is Ca,p 2 Q.'admissible with respect to /x. 

(hi) The bilinear form .A©^^ is closed on L‘^{TL). 


To see that Theorem 2.13 applies to a large class of rough domains, let us recall that fl is said 
to have the kk^’^-extension property if for every u G there exists U G kk^’^(M'^) such that 

U\q, = u a.e. In that case, by [Ml Theorem 5], there exists a bounded linear extension operator 
E from iT^’^(n) into Vk^’^(M'^). It is worth emphasizing that for extension domains the spaces 
W^’^(n) and kk^’^(n) coincide. The following important result is proven in [591 Lemma 4.7]. 


Theorem 2.14. IfEldMJ^ (N > 2) has the -extension property, then dEL is Cap 2 ,n'admissible 
with respect to y = 'H'I'qq, the d-dimensional Hausdorff measure of dEl, for any d G {N — 2, N). 





REACTION-DIFFUSION EQUATIONS 


13 


Furthermore, if /x is an upper d -Ahlfors measure on dfl for some d G {N — 2, N) n (0, N), then by 
[T3l Remark 2.20], dQ is Cap 2 j^-admissible with respect to /x. In particular, since bounded Lipschitz 
dom ains and, for example, the Koch snowflake and the ramified tree domains (see Examples |2 . 3| and 
respectively), have the VF^’^-extension property then their boundaries are Cap 2 ^-admissible 


2.4 


with respect to a. For the Koch snowflake and the tree domain, dQ is also Cap 2 ^-admissible with 
respect to where d = ln(4)/ ln(3) for the snowflake and d = — ln(2)/ ln(a) for the tree domain, 
respectively. Some examples of open sets, whose boundaries are not Cap 2 ^-admissible with respect 
to a are contained in [ 8 l Examples 1.5, 1.6, 4.2, 4.3]. 

2.3. The Laplacian with nonlocal Robin boundary conditions. In this subsection we define 
a realization of the Laplace operator with nonlocal Robin boundary conditions on L^(n). We also 
deduce some properties of this operator that are needed to prove our main results. 


Throughout the remainder of this article, we assume that II C is a bounded domain with 
boundary dQ, /x is a regular Borel measure on dQ such that its domain F^ defined in (2.2) is 


Cap 2 Q-admissible with respect to /x. Then, by Theorem 2.13, the bilinear form defined in 

( |2.11 ) is closed on L^(II). Let be the closed, linear, self-adjoint operator associated with 

in the sense that 

f-D(Ae,^) := {u G Ve’^(II,(9II,/x), 3v e Ae,,,{u,ip) = V g Ve’^(II, dH, /x)} 

[Ae^^xx :=v. 

( 2 . 12 ) 

The operator can be described explicitly as follows: 

D{Aq^^) = {xx G V 0 ^(II, 911, /x) : Axx G T^(II), dyuda + {u + &fj,{u))d^ = 0 on 911}, 

and 

Ae,^iU = -Au. 

We note that if /x = a or, more generally, if fi is absolutely continuous with respect to a (in 
particular, dfj,^ = b{x)dax, for some nonnegative bounded measurable function b{x) G L°° (9II)), 
then 


D { A0 ^„) = {u € V 0 ^(II, 9II, a ) : Au G L^(II), dyU + b{x)u + a{x, y)Qa{u) = 0 on 9II}, 

which corresponds to the classical (or usual) nonlocal Robin boundary conditions. Here a {x, y) = 
b (x) b (y) is a bounded measurable function on dQ x 9II such that dy^d^y = a (x, y) daxday. We also 
notice that the description (2.9) for V 0 ^(n, 911, /x) shows that one always has a Dirichlet boundary 
conditions on the part F^. Clearly the operator defined above, satisfies 


Ae^y : vlf{n,dn,y) ^ (V4’"(II,9II,/x))*, 


1 , 2 / 


(2.13) 


that is, AQ^y maps V 0 ^(II, 911,/x) into (V 0 ^(II, 911,/x)) 


1 , 2 / 


Throughout the remainder of this article, if /x is absolutely continuous with respect to cr, we will 
simply say that y = a. We have the following result. 

Theorem 2.15. Let F^ he n'admissible with respect to /x. Assume that either y = a or that 
V 0 ^(II, 911, /x) ^ L^'^(II) for some q = q {N, II) > 1 if y A (2-14) 

Then the following assertions hold. 
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(a) The operator generates a strongly continuous semigroup of contractions 

on L^( 0 ) which can be extended to contraction strongly continuous semigroups on for 

every p G [ 1 , oo). 

(b) The operator Aq^^ has a compact resolvent, and hence has a discrete spectrum. The spec¬ 
trum of Aq^^ is an increasing sequence of real numbers 0 < Ai < A 2 < • • • < A„ < ..., 
that converges to +cx). Moreover, if Un is an eigenfunction associated with Xn, then Un G 
D{Ae^^)nL^{n). 

(c) Assume also that either Cap 2 ^f 2 (r“) > 0 or > 0. Then 0 G p{Aq^^), the resolvent set 

of Aq^^, and there exists a constant C > 0 such that, for every u G Vq'^ dQ, p), 

\\u\\l^^ <CAe,i,{u,u). (2.15) 

(d) Denoting the generator of the semigroup on U’{Q) by Ap^Q^^^, so that = ^ 2 , 0 ,/^? then 
the spectrum of Ap^Q^^ is independent of p. The spectral bound 

s = inf {Re A : A G a 

of Ap^Q^fj, is an algebraically simple eigenvalue. It is the only eigenvalue having a positive 
eigenfunction. 

(e) For each 0 G (0,1], there holds D{Aq C (0) provided that 0 > 7 . Here, 7 = N/2 if 
p = a and 7 = q/2{q — 1) if p A 


Proof, (a) This part follows from |59l Corollary 5.3] (where the assumption (2.14) is not needed). 

(b) First, we notice that since we have assumed that F^ is Cap 2 j^-admissible with respect to 
/i, we have from Theorem 2.13 that the continuous embedding V 0 ^(n, 911,/i) ^ L^(ll) is also an 
injection. Moreover it follows from (2.6) if p, = a and the assumption ( |2.14 ) if 7 ^ cj, that the 
embedding V 0 ^( 11 , 911,/r) Lp‘{Tl) is compact since H is bounded. Since ^ 0 ,^ is a nonnegative 
self-adjoint operator and has a compact resolvent, then it has a discrete spectrum which is an 
increasing sequence of real numbers 0 < Ai < A 2 < • • • < A„, •.., that converges to -|-oo. 

Now, let Un G V0^(11,911, ^) be an eigenfunction associated with Xn. Then, for every v G 
V^’2(11,911,/x), 


VMri'Vudx+ / UnVdp 
! iTa 


+ 


Ks (x, y) {Un{x) - Un{y)){v{x) - v{y))dp^dp =Xn UnVdx. 

Jn 


This equality means that A^^pUn = XnUn- Let a > 0 be a real number. Since a G p{—Aq^p), we 
have that al + is invertible. From A^^pUn = XnUn we have that 

— (ciL T Aq p') {Xn T Ot)Un — (A^ T Ol){cxI T Aq p^ (Unf 

By |59l Theorem 5.4], the semigroup (e“*^®’'^)t>o is ultracontrative in the sense that it maps L^(H) 
into L°°(ll) (this follows from ( |2.6| ) if ;U = a, and from the assumption (2.14) if p A ^)- More 
precisely, there is a constant C > 0 and 7 G (0,1) such that for every / G L^(ll) and t G (0,1), 


l|e-“'’"7lloo,n < Ct-' ||/||j„ (2.16) 

where 7 = N/2 if p = a and 7 = q/2{q — 1) if p A cr. Since for every / G L^(H) and a > 0, 

roo 

{al + Aq^p) ^f = 

Jo 
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it follows from (2.16) that there exists a constant M > 0 such that 

< M{\n + a)\\Un\\2,n- 


Un 


This completes the proof of part (b). 

(c) Assume to the contrary that 0 ^ p{AQ^y) (i.e., 0 is an eigenvalue for Then there 

exists a nonzero function u G D{AQ^fj) such that Aq^^u = 0, that is, As,^iiu,v) = 0, for every 
V G V 0 (n,(9n,/i). In particular, we have 


Ae^fj,{u,u) := / \'Vu\‘^dx + / \u\‘^dfjL + 


Ks{x, y)\u{x) - u{y)\^dp.^dp = 0 . 


This shows that |Vn| = 0. Since is bounded and connected, we get that u = Ci for some constant 
Cl- Since u = 0 q.e. on r“ (as u belongs to V 0 ^(fl, cAI,/i)), there are two alternatives. First, 
CaP 2 ,o(rjf) > 0 in which case u = Ci ^ D{Aq^^), whence a contradiction. On the other hand, from 
the facts that u = Ci, we get from the previous equation that C'i/i(r^) = 0. The second alternative 
/i(r^) > 0 then implies that u = Ci = 0, again a contradiction. We have shown that 0 G p{AQ^y). 
Now the estimate (2.15) follows from standard properties of semigroups (see, e.g., naEO]) , since 
we have just shown that there exists a constant w > 0 such that \\e~^^^'>^\\c(L'^{n),L'^{n)) < for 
every t > 0 . This completes the proof of part (c). 

(d) Let be the generator of the semigroup on LP{Q). Since Aq^^ = ^ 2 , 0 ,^ has a compact 

resolvent and is bounded, it follows from the ultracontractivity that each semigroup has a compact 
resolvent on L’p[Q) for p G [l,oo]. Now it follows from [23L Corollary 1.6.2] that the spectrum of 

is independent of p. The last assertion in (d) follows from |18l Proposition 5.9 and Corollary 

5.10]’. 

(e) Since I + A©^p is invertible we have that the (flj-norm of (/ + Aq^^)^ defines an equivalent 
norm on D[Aq and 


(/ + Ae^p) ^ — 


r(0) 




fdt. 


Using (2.16) for t G (0,1) and the contractivity of e for t > 1, for u G ^^(Aq ^), we deduce 


I^IIl°°(o) — ^ 11^ 


^ '^dt + C \\u 


I,m) 


fC 


-t 


dt. 


The first integral is finite if and only if 7 < 0. The proof of theorem is finished. 

Remark 2.16. 


□ 


2.15. is actu- 


(a) We remark that either Cap 2 ^o(r“) > 0 or piT^) > 0, in part (c) of Theorem 
ally not an assumption about the regularity of the open set Q, but it is only an assumption 
about the measure p. Roughly speaking, it means that the trivial ehoice p = 0 (the zero 
measure) is not allowed. 

(b) In the case of domains with Lipschitz continuous boundary dH, further interesting spectral 
properties for the operator Aq^^ can be also derived (see [33] ). 

(c) Assuming that either Cap 2 ^Q(r]]°) > 0 or p{T^) > 0 is satisfied, we deduce from (2.15), 

1/2 




v^'^{n,dn,u) 


2,0 + ll^llL2(r^) +Afs,iiiu)^ 


(2.17) 


defines an equivalent norm on the space V 0 ^(U, cAI,/i). In particular, under the same as¬ 
sumption, we also have that 

1/2 

■“ (ll'^'*^ll 2 ,o + (2.18) 
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defines an equivalent norm on the space h)- 

(d) In the case when 9 G (0,1] and pL = a, the embedding D{Aq C (Q) holds provided 
that N < 26. T his consideration is optimal for an arbitrary open subset Q C in view 
of Remark 2.1[\ Finally, when has the -extension property of Sobolev functions, we 
have q = N/ (N — 2) for N > 3 and q G (1, oo) for N = 2. In that ease D{Aq C (H) 
provided that N < AO if N > 3. In a smooth situation (say when 511 G C'^), one has 
D{AQ^fj) C (H) and thus D{Aq C It is well-known that the embedding 

(H) C (H) holds in the smooth setting provided that N < A6. This shows that the 
stated embedding D{Aq ) C (H) appears to be optimal in the non-smooth setting in this 


sense. 


It is important to observe that the operat or depends on the choice of the measure p, on 
5H. In order to see this, consider Example 2.3 when H is the Koch snowflake in and take 
/I = (T = For this choice, a is locally infinite on 511 so that there are no Robin boundary 

conditions for H, on account of the description (2.8)-(2.9). Indeed, we have [ 8 l Theorem 2.3], 

R (-^ 0 , 0 -) = {u G (H) : u = 0 q.e. on 5H} = W^''^ (H) 

and, consequently, ^e^o- coincides with the realization of the Laplace operator —A with homoge¬ 
neous Dirichlet boundary conditions on 511. However, with a different choice p = (where 

d = In (4) / In (3) is the Hausdorff dimension of 511), we have the generalized Robin boundary con¬ 
dition ( 1 . 2 ) on 5H since r“ 


= 0 (see Example 2.3). This is also the case for the tree-domain H 


of Exarnple 2.4 in which case Aq^^ coincides exactly with the realization of the Laplace operator 
with homogeneous Dirichlet boundary conditions on Tj (the small structures), and nonlocal Robin 
boundary conditions on 511 \Tf. However, recall that the fractal boundary Tj of 11 is a d-set, with 
d = — In (2) / In (o). Thus, the natural measure p of 511 should be once again the d-dime nsion al 
Hausdorff measure p = 'H'I'qq so that the generalized (nonlocal) Robin boundary condition (1.2) is 
still realized on the portion Tj. In fact, a more appropriate choice for p is 

O', 


h = 


on Tf, 
on 511 \ T f, 


so that one has the usual (nonlocal) Robin boundary conditions on the portion 511 \ Tf, and the 
generalized (nonlocal) Robin boundary conditions onTj. In the general case when H C M.^, such 
arguments also hold if 511 is a d-set with d G {N — 1, A^]. 

Finally, it is worth emphasizing that the assumption ( |2.14 ) in Theorem 2.15 is satisfied with 
q = N/{N — 2) if, for example, 11 has the VE^’^-extension property. On the other hand, if H C 
is a non-Lipschitz bounded domain whose boundary admits a finite numbers of outward or inward 
Holder cusps points and exponent at cusp points 0 < 7 < 1, then H does not enjoy the extension 
property. However, for such domains H we still have the following continuous embeddings: 


V0’^(ll,511,/i) -A -A with q := ^j ^ > 1 , 


(2.19) 


see [281129j . Many examples of domains with cusps, which are the simplest non-Lipschitz domains 
used in the applications, can be found in [13 Chapters 7-8]; see, for instance, Example 2.17 below. 

that 


we have from Definition 


2.1 


For Holder C''’'-domains and for the measure p = a = 
r“ = 0, whereas = D (a) = 511. Thus, it is clear that in such cases, the usual (nonlocal) Robin 
boundary condition ( 1 . 2 ) is satished on To- = 511. 
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Figure 3. Cusp-shaped domain 


Example 2.17. (see Figure 3) Let Like C N >2 he the following set: 

Fthc := ja; = (xi,X2, ...,xn) G : xat G (0, L) , 0 < > 


for some L,l > 0 and 7 G (0,1). The inequality (2.11^ is fulfilled for Ll = and N >2. Generally, 

the Sobolev inequality is satisfied in any Holder domain Cl of elass C'^, 7 G (0,1). By a 

(bounded) Holder C'^-domain we mean the following: there exist a finite number of balls B {xi,ri), 
i = 1, whose union eontains fl, and for each i = I, there exists a Holder-continuous 

(C'^, 7 G (0, 1)) function Fi : —)• M such that for some coordinate system, B {xi,ri) n is 

equal to the intersection of B {xi,ri) with the region above the graph of Fi. 

3. Well-posedness and regularity 

In this section H C is an arbitrary bounded domain with boundary dCl, is a regular Borel 
measure on dCl such that pL is either the {N — l)-dimensional Hausdorff meausure or another 

’’arbitrary” regular Borel measure, according to the assumptions of Section]^ (see (H^) below). We 
begin this section by stating all the hypotheses on / we need, even though not all of them must be 
satisfied at the same time. 

(HI) / G Cjqp (M) satisfies 

f{s) 


lim inf 

|s|^+oo 


> -A* 


for some constant A* G [ 0 ,( 7 ^), where Cq = C {N,Cl) > 0 is the best Sobolev/Poincare 
constant into (2.5) and (2.14), respectively. 

(H2) / G (M) satisfies 

C/ |s|^ — c/ < / (s) s < C/ |s|^ -|- c/, for all s G M, 

for some appropriate positive constants and some p > 1. 

(H3) / G Cl (M) satisfies 

/ (s) > —C/, for all s G M, 
for some positive constant Cf. 
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Finally, throughout the remainder of this article our assumption about the measure /U will be as 
follows. 

(H;.) Assume that F^ is Cap 2 ^Q-admissible with respect to /i, and either Cap 2 ^o (r“) > 0 or 

/i(F^) > 0. 


In what follows we shall use classical (linear/nonlinear semigroup) dehnitions of strong and 
generalized (weak) solutions to (1.1)-(1.3). “Strong” solutions are defined via nonlinear semigroup 
theory for sufficiently smooth initial data and satisfy the differential equations almost everywhere 
in f > 0. “Generalized” or weak solutions are defined as (strong) limits of strong solutions. We 
first introduce the rigorous notion of (global) weak solutions to the system (!.!)-(1.3). Throughout 
the remainder of this article the solution of the system is a function that depends on 

the variables t G [0, oo) and x G but in our notation we sometime omit the dependence in x. 

Definition 3.1. Let uq G L?‘{LL) be given and assume (H2) holds for some p > 1. The function 
u is said to he a weak solution of EH - ( ff .gp if, for a.e. t G (0,r), for any T > 0, the following 
properties are valid: 

• Regularity: 

f n G (0, T; (Q)) n RP ((0, T) x Ll) n L‘^{0, T; Ve’^(0, dQ, p)), 

[ dtu€ L^o,T-(y^/{n,dn,p))*) + LP ((o,r) x n), 


(3.1) 


where p = p/ {p — 1 ). 

Variational identity: for the weak solutions the following equality 


(3.2) 


{dtu (t ), () A 0 ,f,{u {t) ,f) + {f {u (t)) ,0 = 0 

holds for all ^ G V 0 ^(n, (911,//) n RP {Q), a.e. t G (0,r). Finally, we have, in the space 
L?' (n), M ( 0 ) = uq almost everywhere. 

Energy identity: weak solutions satisfy the following identity 


1 „ 

2 11 ^ 


In + 


V 1 

{u (s), u {s)) (/ {u (s)), u (s))] ds = - ||n 


/o 


2 

2,0 


(3.3) 


Remark 3.2. Note that by (3.1), u G C^ieak ([0, T] ; L^ (11)). Therefore the initial value u (0) = uq 
is meaningful when uq G L^ (H) . 

Finally, our notion of (global) strong solution is as follows. 

Definition 3.3. Let uq G L°°(11) be given. A weak solution u is ’’strong” if, in addition, it fulfills 
the regularity properties: 

^ e 0“((o, T]-,L\n)) n C ([0, T] ; 

such that u (t) G D{Aq^^), a.e. t G (0, T), for any T > 0. 


m), 


(3.4) 


This section consists of three main parts. At first we will establish the existence and uniqueness 
of a strong solution on a hnite time interval using the theory of monotone operators, and exploiting 
a (modified) Alikakos-Moser iteration-type argument to show that the strong solution is actually a 
global solution. In the second part, we will show the existence of (globally-defined) weak solutions 
which satisfy the energy identity (3.3) and the variational form (3.2). We shall also discuss their 
uniqueness. Finally, by using the energy method combined with a refined iteration scheme, we show 
that the weak solution is actually more regular on intervals of the form [d, oo), for every d > 0 . 
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3.1. Strong solutions. Now we state the main theorem of this (sub)section. 

Theorem 3.4. Assume that the nonlinearity f obeys (HI) and /i satisfies (H^). For every uq G 
L°° (n), there exists a unique strong solution of (1.1)-(1.3) in the sense of Definition 3.3. 

Proof. Step 1 (Local existence). Let uq G L°° {Ft) C (fl) = D(Aq^p() ^ ^ = D{(pQ^^) ^ \ By 

Theorem 2.131 is proper, convex and lower semi-continuous on the Hilbert space LP‘{FL) (cf. 

we know that (= —dipQ^^) generates 


also [59 


heorem 5.2]). Moreover, from Theorem 


2.15 


a strongly continuous (linear) semigroup (e of contraction operators on L‘^{Fl). Finally, 

(Appendix), is non-expansive on (H), 


by Theorem 


6.1 


|e < ||iio|loo,D > t > 0 and no G (11), 


(3.5) 


and Aq^^ = is strongly accretive by Theorem 2.15-(c), see (2.15). Thus, we can give an 


operator theoretic version of the original problem (1.1 )- p.3p ; 

dtu = -Ae^^u - / (n), (3.6) 

with / as a locally Lipschitz perturbation of a monotone operator, on account of the fact that 
/ G Cjqp (M). We construct the (locally-defined) strong solution by a fixed point argument. To this 
end, fix 0 < T* < T, consider the space 

= {n G C([0,rl ;L~ (H)) : \\u{t)\\^^^ < il*} 

and define the following mapping 

S (n) (t) = e-*^®’'^no - / (n (s)) ds, t G [0, T*] . (3.7) 

Jo 

Note that when endowed with the norm of C ([0, T*] ; (11)), is a closed subset of 

C* ([0, T*] ; (h)), and since / is continuously differentiable, S (u) (f) is continuous on [0,T*]. 

We will show that, by properly choosing T*,R* > 0, S : Xt*,r* —>• Xt*,r* is a contraction mapping 
with respect to the metric induced by the norm of C ([0, T*] ; L°° (H)) . The appropriate choices for 
T*,R* > 0 will be specified below. Firs t, w e show that u G Xt*,r* implies that S (rt) G Xt*^r*, 
that is, S maps Xt*,r* to itself. From (3.5) and the fact that / G (M), we observe that the 
mapping S satisfies the following estimate 


S(u(())IL,n<IKIUn+ / <!-<‘-l-’<».^{/(0) + {/(u(s))-/(0))) 

JO 

<\\ML,n + t{\fi0)\+Qf {R*)R' 


oo,r 2 


ds 


for some positive continuous function Qji which depends only on the size of the nonlinearity f'. 
Thus, provided that we set R* > 2 ||wo||oo q, we can find a sufficiently small time T* > 0 such that 


2 T 


(|/(0)|+Q/ {R*)R*)<R*, 


(3.8) 


in which case S {u{t)) G Xt*,r*, for any u{t) G Xt*,r*. Next, we show that by possibly choosing 
T* > 0 smaller, E : Xt*,r* —>■ Xt*,r* is also a contraction. Indeed, for any u\,U2 G Xt*,r*, 
exploiting again (3.5), we estimate 

rt 


|E(i.i(i))-E(.«(i))IUn<Q;.(R*) f 

< tQ(R*) ||ui - . 


oo,n 


ds 


(3.9) 
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This shows that S is a contraction on Xt*,r* (compare with (3.8)) provided that 
2T*{\fm+Qf,{R*)R*) _2T*{\fm+Qf,{R*))R* 


max 




-,T*Qf,{R*)\ = 


R* 


< 1 . 


Therefore, owing to the contraction mapping principle, we conclude that problem (3.6) has a unique 
local solution u G Xt* ,_r* . This solution can certainly be (uniquely) extended on a right maximal 
time interval [0,Tmax)) with T^ax > 0 depending on HuollcoO) such that, either Tmax = oo or 
Tmax < oo, in which case lim^/'Tmax 11^ (Olloo n ~ Indeed, if Tma x < oo and the latter condition 


does not hold, we can find a sequence tn T^a x such that \\u (tn)|loo n — This would allow us to 
extend u as a solution to Equation (3.6) to an interval [0, tn + 6), for some <5 > 0 independent of n. 
Hence u can be extended beyond Tmax which contradicts the construction of Tmax > 0. To conclude 


G G lTh2 [6, Tmax; (II)) n (<5, Tmax; (H)) • 


that the solution u belongs to the class in Definition 3.3, let us further set G (t) := —/ {u (t)), for 
u G C ([0,Tmax) ;T“ (D)) and notice that u is the ’’generalized” solution of 

dtu + = G{t), tG [ 0 , Tmax), (3.10) 

such that u (0) = uq G L°° (D) C = D{Aq^^). By Theorem 6.2 (Appendix), the ’’generalized” 

solution u has the additional regularity dtu G [S, Tmax ; (H)) , which together with the facts 
that u is continuous on [ 0 , Tmax ) and / G (M), yield 

(3.11) 

Thus, we can apply Theorem |6.3| and Corollary |6.4| (see Appendix), to deduce the following regu¬ 
larity properties for u : 

UGL°° ([<5,Tmax);T(Ae,^)) n lyhoo ([5^rmax);T2 (D)) , (3.12) 

such that the solution u is Lipschitz continuous on [5, Tmax), for every > 0. Thus, we have obtained 
a locally-defined strong solution in the sense of Definition |3.3[ As to the variational equality in 
Definition 3.1, we note that this equality is satisfied even pointwise (in time t G (0, Tmax )) by the 
strong solutions. Our final point is to show that Tmax = oo, because of condition (HI). This ensures 
that the strong solution constructed above is also global. 

Step 2 (Energy estimate) Let m > 1 and consider the function Em ■ (0, oo) —)> [0, oo) defined 
by Em{t) := ||^i(l)|lm+i n- First, notice that Em is well-defined on (0, Tmax ) because u is bounded 
in D X (0, Tmax ) and because Gl has finite measure. Since u is a strong solution on (0, Tmax ), 
see Definition 3.3 (or (3.12)), recall that u is continuous from [0,Tmax) L°°{Gl) and Lipschitz 
continuous on [5, Tmax) for every <5 > 0. Thus, u (as function of t) is differentiable a.e., whence, the 
function Em{t) is also differentiable for a.e. t G (0,Tmax) • 

For strong solutions and t G (0,Tmax), integration by parts procedure yields the following stan¬ 
dard energy identity: 


\ d f 

2 ~dt^^ (*) {u{t),u{t)) + y^/(^ (I))^(I) dx = 0. 

Assumption (HI) implies that 

/ (s) s > - Cf, 

for some Cf > 0 and for all s G M. This inequality allows us to estimate the nonlinear term in 
(3.13). We have (by using the equivalent norm (2.17[)) 


(3.13) 


(3.14) 




< Cf\n\+X,Ei (t), 


(3.15) 
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where |r2| denotes the A^-dimensional Lebesgue measure of fi. In view of (2.17) and Gronwall’s 

(3.16) 


inequality, (3.15) gives the following estimate for t G (0, Tmax), 

rt 


ll« (0112,0 + 2(^0-A*) / h('S)llvP(n,ao,/.)^® - Il'“oll2,oe ^ +<+(/, |11|), 

for some constants p = p {N, Q) > 0, C (/, |11|) > 0. 

Step 3 (The iteration argument). In this step, c > 0 will denote a constant that is independent 
of t, Tjnax, rn, k and initial data, which only depends on the other structural parameters of the 
problem. Such a constant may vary even from line to line. Moreover, we shall denote by Qr {m) a 
monotone nondecreasing function in m of order r, for some nonnegative constant r, independent 
of m. More precisely, Qr (m) ~ cm"^ as m —+oo. We begin by showing that Em (t) satisfies a local 



following inequality: 




VImI 2 


dx + 


\u 


m+1 


dp 


(3.17) 


< Qi{m + l) [Em{t) + {Em{t))'^+^ ^ 
for some constant 7 > 0 independent of m, t and Tmax- Next, set ruk + 1 = 2^, /c G N, and define 

(3.18) 


- max* 
| 2 ^ 


Mfc := sup / \u{t,x)\ dx = sup Em^{t) ■ 

IG(0,Tinax) ^ le(0,Tmax) 


Our goal is to derive a recursive inequality for using (3.17). In order to do so, we define 

1 


Pk ■= 


kklk 'kTlk—l 


q {1 + mk) - {1 + nik-i) 2q - 1 


< 1 , Qk-='^-Pk = ‘^ 


q-l 

2 q-l 


since q > 1. We aim to estimate the terms on the right-hand side of (3.17) in terms of the 
(fl)-norm of u. First, Holder inequality and Maz’ya inequality (i.e., VF 2 2 (H, clH,/r) C 
(H), for some q = q{N,Q) > 1, see (2.6) with q = N/ {N — 1) and p = a, and (2.14) when 

(3.19) 


p ^ (j) yield (by using the equivalent norm (2.18)) 


u 


in 


u 


'n 


< c 


in 




V 2 


Pk 


u 




) Qk 




dx + 


\u 


l+rrik 


dp 


\u 


l+nik-i 


dx 


Qk 


with Sfc = p^q = q/ {2q — 1) G (0,1). Applying Young’s inequality on the right-hand side of (3.19), 
we get for every e > 0 , 


Qi {mk + 1) / 
Jn 


u 


l+ruk 


dx < e 


V \u\ 




dx + £ \u 
JTu. 


+Qa {pnk + ^^ 


u 




dp 

max{ 2 , 2 fc} 


(3.20) 


dx 


for some a > 0 independent of k, where in fact Zk '■= qk/ {^ ~ Sk) = 2 , for each k, since mk -|-1 = 2 ^. 
In order to estimate the last term on the right-hand side of (3.17), we define two decreasing and 





























22 


CIPRIAN G. GAL AND MAHAMADI WARMA 


increasing sequences and , respectively, such that 

Qk'mk 


mk + 1 , 

Ik = ^- and Wk ■= 




nik (1 - Sfc) + 1' 


and observe that they satisfy 


1\ 2(g-l) 

1 < 4 < 2 2 - - , <Wk<2 


qj - 2 


for all A: G N (in particular, tcfc —>■ 2 as A —)• oo). The application of Young’s inequality in (3.20) 
yields again 


m-k 

Qi (wifc + 1) 


< e 


l + rrik 

V 


dx + 


\u 


l+ruk 


d/i (3.21) 






dx 


Wk 


for every e > 0 , where now 


(mfc + 1)^ {l + rukY 


with Pk' '■= i'krik + 1) / {ni k (1 — Sk) + 1) —^ {2q — 1) / {q — 1), as fc —>■ oo. Hence, inserting (3.20), 
(3.21) into equation (^3.17|, choosing a sufficiently small 0 < e < eo := 3 Riia (7)1)) and simplifying, 
we obtain for t G (0,. 


^ J |u (t, dx 

< Qs {mk + 1) |r/ 


|^+™'=-idx 


1 +m^ 

V 

2 


dx + 


\u 


mk+l 


d/i 


(3.22) 


for some positive constant <5 > 0 independent of k. 

1 +mfc 

Next, we have |u| 2 g H 2^2 {Q,dQ,^) (recall that u{t) G VF 2^2 (12)<9n,/i) n L°° (Q) , for a.e. 

t G (0,Tmax))) so that by using the equivalent norm ( 2.18| ), we can find another positive constant 
c = c (n, N), depending on Q and N, such that 


l + rrik 

V 


dx + 


irrife+l 


d/Li > c / \u 

Jn 


rrik + l 


dx. 


We can now combine (3.23) with (3.22) to deduce 
d 


~dt + ^ j^\u{t,x)\^*‘ dx <Q 5 


(3.23) 


(3.24) 


for t G (OjTmax) • Integrating (3.24) over (0, t), we infer from Gronwall-Bernoulli’s inequality |17L 
Lemma 1.2.4] that there exists yet another constant c > 0, independent of k, such that 


Mfc < max i I \uq\^ dx, )■ , for all A: > 2 


(3.25) 

00,0) — 1’ 


On the other hand, let us observe that there exists a positive constant Coo = Coo(||uo 
independent of k, such that ||uo|l 2 feo ^ ^oo- Taking the 2^-th root on both sides of (3.25), and 
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defining Xk := sup^g(Q ||n (i)|l 2 *: n > easily arrive at 

Xk < max I Coo, Xfc_i| , for all k >2. 


(3.26) 


By straightforward induction in (3.26) (see [71 Lemma 3.2]; cf. also [T71 Lemma 9.3.1]), we finally 
obtain the estimate 


sup 

te(0,Tma 


0 


I^WIIoo.D < 


lim Xk < c max 
k^+oo 


a 


oo, sup 

t£{0,Tm3 


0 


|n(t)|| 


2,0 


(3.27) 


It remains to notice that (3.27) together with the bound (3.16) shows that jju q is uniformly 
bounded for all times t > 0 with a bound, independent of Tmav depending only on lluoUooO) the 
’’size” of the domain and the non-linear function /. This gives T^ax = +oo so that strong solutions 
are in fact global. This completes the proof of the theorem. □ 


Remark 3.5. Strong solutions to Eqs. (1.1)-(1.3) exhibit an improved regularity in time, we have 

u€C ([0, T ]; (11)) n C{{6, OQ, /x)), (3.28) 

for any T > 6 > 0. This follows from the fact that the nonlinear function / is continuously 


differentiable. Note that the second regularity in (3.28) is a consequence of the first one, the time 


regularity in (3.4) (see Definition 3.3) and the variational identity (|3.2p . 
The following result is immediate. 


Corollary 3.6. Let the assumptions of Theorem 3-4 be satisfied. The reaction-diffusion system 
(1.1)-(1.3) defines a (nonlinear) continuous semigroup T (t) : L°° (D) —)• (D), given by 

T (t) uq = u (t ), (3.29) 


where u is the (unique) strong solution in the sense of Definition 3.3. 

3.2. Weak Solutions. We aim to prove the existence of weak solutions in the sense of Definition 


3.1 For this. Theorem |3.4| proves essential in the sense that we can proceed by an approximation 
argument, which we briefly describe below. We consider, for each e > 0, the following system 

dtUe — Xue + / (Uf,) = 0, in D X (0, oo), 

duUeda -|- (tie -|- 0 ^ {ufij)dp = 0 , on dD x ( 0 , oo), 

subject to the initial condition 

Ue (0) = UOe, 

where uoe £ L°° (D) n V 0 ^( 11 , dD, ff) such that 

UOe Uq = u ( 0 ) in (D). 


(3.30) 


(3.31) 


(3.32) 


Then, by Theorem 3.4, under the assumption that / satisfies (HI), the approximate problem 


(3.30)-(3.31) admits a unique strong solution with 

u, E ((0, T4-,L\n)) n C ([0, T,] ; (D)) n Lg’, ((0, T,]; D (H©,^)), 


(3.33) 


for some T* > 0 and each e > 0, such that both (3.2) and 


are satisfied even pointwise in time 
t E (0, T*). The advantage of this construction is that now every weak solution can be approximated 
by the strong ones and the rigorous justification of our estimates for such solutions is immediate. 


We shall now deduce the first result concerning the solvability of problem (1.1)-(1.3) with the 


new assumption (H2). Note that a function that satisfies (H2) also satisfies the assumption (HI). 
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Theorem 3.7. Assume that the nonlinearity f obeys (H2) and fj, satisfies (H^). Then, for any 
initial data uq £ , there exists at least one (globally-defined) weak solution 

u£C{[0,T]-,L‘^{n)) 

in the sense of Definition\3.1\ 


Proof. We divide the proof into two main steps. For practical purposes, C will denote a positive 
constant that is independent of time, T*, e > 0 and initial data, but which only depends on the 
other structural parameters. Such a constant may vary even from line to line. 

Step 1. (The main energy estimate). Since / obeys (H2), then (HI) is satished and by Theorem 


3.4, the approximate problem (3.30)-(3.31) has a strong solution that also satisfies the weak for¬ 


mulation (3.2). Thus, in view of (3.33) the key choice {t) in (3.2) is justified. We have the 

following energy identity 

\ d f 

2 )fi\\'^'^(*)\\ 2 ,n + -^e,ti{ueit),u,{t))+ / f {ueit))u, (t) dx = 0, 


for all t £ (0,T*). Invoking assumption (H2), we infer 

|2 


11“^ (^)ll|o + \\Ue 




;,n<c\n\ 


(3.34) 


where we have used the equivalent norm (2.17). We can now integrate this inequality over (0,Ti,) 
to deduce 

\K{t)\\l,n+ [ {CQ\\ue{s)\\li^2^^g^^^^+Cf\\ueis)fp^^'jds<\\ue{0)\\l^^e-P^ + C (3.35) 
^ 0 

for all t £ (0,T,), for some p > 0 independent of e > 0. Incidently, by (3.32) this uniform 
estimate also shows that we can take T* = oo, i.e., the weak solution that we construct is in fact 
globally-defined. On account of ( 3.35[ ), we deduce the following uniform (in e > 0) bounds 

u, £ (0, T; (H)) n (0, T; V (H)), (3.36) 

u, £ L2(0,r;Ve’^(H,(9O,/i)), 

for any T > 0. Hence, by virtue of (2.13) and ( |3.36 ), we also get 

£ l2(0,T;(V^’'(H,5H,//))*), f {uf) £ L^' (H)), (3.37) 

uniformly in e > 0. Here, recall that Aq^^ is the positive self-adjoint operator associated with the 
bilinear form .40,^ (see Section 2.2). By (2.13), we can interpret (3.2) as the equality 

- dtUe = Ae^^Ue -A f {Ue) (3.38) 


in L^(0, T; (V 0 ^(H, (9H,/i))*)-|-(0, T; Z/P (H)), where p =p/{p—l) (see also (3.6)). 

Step 2. (Passage to limit). From the above properties (3.36)-(3.37), we see that there exists a 
subsequence (sHU denoted by {ue}), such that as e —)• O"'', 

Ue ^ u weakly star in L°° (O, T; (H)) , 

u weakly in PP ( 0 , T\LP {Pl)), ^ ^ ( 3 . 39 ) 

dtUe —)• dtu weakly in L^(0, T; (V 0 ^(H, 50,/r))*) + 17 {0,T;LP (O)). 

Since the continuous embedding V 0 ^(O, 50, fi) (0) is compact, then we can exploit standard 

embedding results for vector valued functions (see, e.g., m). to deduce 

u strongly in (O, T; (O)) . 




(3.40) 
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By refining in (3.40), converges to u a.e. in VL x (0,T). Then, by means of known results in 
measure theory m, the continuity of / and the convergence of ( |3.40[ ) imply that / (n^) converges 

/ I I I I '-' 

weakly to f (u) in ((0,r) x 14), while from (3.36)-(3.37) and the linearity of we further 

see that 


Ae,^iUe Ae^fj,u weakly in L^(0, T; (V 0 ^(n, 514, ^))*). (3.41) 

We can now pass to the limit as e —)• 0 in equation ( |3.38 ) to deduce the desired weak solution u, 
satisfying the variational identity (3.2) and the regularity properties (3.1). In order to show the 


energy identity (3.3), we can proceed in a standard way as in |161 Theorem II. 1.8], by observing that 

any distributional derivative dtu {t) from V'{[Q, T ]; (V0^(14, 514, ^))* + (H)) can be represented 

as dtu (t) = Zi (t) + Z 2 (t ), where 

(t) := -Ae,^u (t) G L\0, T; (V^’2(14,514, /x))*), ^2 (t) := -/ (u (t)) G (0, T; L^' (14)), 

which are precisely the dual of the spaces L^(0, T; (V0^(14, 514,/x))), and LP(0, T; (14)), respec¬ 
tively. In particular, we obtain that every weak solution u £ C ([0, T] ; (H)) , and that the map 

1 1 —>■ ||u (t )||2 Q is absolutely continuous on [0, T], such that u satisfies the energy identity 


\\u (i)ll 2 ,n + {Ae,^,u (t), u (t)) + (/ {u (t )), u (t)) = 0 , 


(3.42) 

□ 


a.e. t G [0, T], whence (3.3) follows. The proof of Theorem 3.7 is finished. 

As in the classical case [m ED |53], we can prove the following stability result for the class of 


weak solutions constructed by means of Definition 3.1 


Proposition 3.8. Let the assumptions of Theorem 3.1 he satisfied, and in addition, assume that 
(H3) holds. Then, there exists a unique weak solution to problem (1.1)-(1.3), which depends con¬ 
tinuously on the initial data in a Lipschitz way. 

Proof. As usual, consider any two weak solutions ui,U 2 , and set u{t) = xxi (t) — U 2 (t). Then, 
according to (3.3) (cf. also (3.42)), u{t) satisfies the identity 

\\uit)\\ln+Ae,^, {u {t),u{t)) = - {f iui (t)) - / (u 2 {t)),u{t )), 
a.e. t G [0, T] . Assumption (H3) implies (by using the equivalent norm (|2.17[)) 


d „ 
dt 


2 ,n 


+ C\\u 




This yields the desired result by application of Gronwall’s inequality. 


(3.43) 

□ 


As an immediate consequence of this stability result, problem (!.!)-(1.3) defines a dynamical 
system in the classical sense [inillSlEDES]. 


Cor ol lary 3.9. Let the assumptions of Proposition be satisfied. The reaction-diffusion system 
defines a (nonlinear) continuous semigroup S (t) : L^ (14) —>• L^ (14), given by 


S (t) uo = u (t ), 


(3.44) 


where u is the (unique) weak solution in the sense of Definition 3.1 
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3.3. Regularity of solutions. The next proposition is a direct consequence of estimate (3.24) of 
Theorem 3.4, Strong solutions guaranteed by Theorem 3.4 provide sufficient regularity to justify 
all the calculations performed in the proof of Proposition 3.10 below. In this case, at the very end 
we pass to the limit and obtain the estimate even for the generalized solutions. 


Proposition 3.10. Let the assumptions of Theorem 3-4 or Theorem 3.1 be satisfied. Letr > r > 0 
and fix p := T —T. There exists a positive constant C = C (p) ~ (for some i] > 0), independent 
oft and the initial data, such that 


sup ||u 

t>r' 


n ^ <^sup||u(s)||2_f5. 


S>T 


Proof. The argument leading to (3.45) is analogous to [311 Theorem 2.3] (cf. also 


(3.45) 

Theorem 


2.3]). It is based on the following recursive inequality for (t), which is a consequence of (3.24) 
and ( [3T^ -( [33^ : 

sup Em,, {t) < C ( 2 ^) ( sup («) ) , for all fc > 1, (3.46) 

^ ' \S>t^: J 

where the sequence {tfclfcgN defined recursively t^ = tfc-i — ijl/ 2^, k > 1, t^ = t' . Here we recall 
that C = C (p) > 0, h > 0 are independent of k. For the sake of completeness, we report a sketch 
of the argument for ( |3.46 ). To this end, let ((s) be a positive function ( : M+ —>■ [0,1] such that 
C (s) = 0 for s G [O, t — ^/2^] , C (s) = 1 if s G [t, +oo) and \dC/ds\ < 2^/fi, if s G (t — /i/2^, t). We 
dehne (s) = ( (s) Em,, (s) and notice that 

(3.47) 


(s) < ( (s) (s) + ^ -E'mfc (s) 


-C{s) ^^Em, {S) + Ql ( 2 ^)^ 


\u 


l+rrife 


dx. 


The last integral in (3.47) can be estimated as in (3.20) and (3.23). Combining the above estimates 
and the fact that Zj, < Em,,, we deduce the following inequality: 


^Zk{s) + C2’^Zk{s)<C 

as 


( 2 ^) ( sup (s))^, for all s G t- p/2^, + 00 ) . 


(3.48) 


Note that C = C (p.) ~ ^ ^ as /r —)• 0, and C (p) is bounded if /x is bounded away from zero. 
Integrating ( |3.48 ) with respect to s from t — /x/2^ to t, and taking into account the fact that 
Zk {t - Ai/2^) = 0, we obtain that Em,, {t) = Zk (t) < C {2^y (sup^>i_^/ 2 fc Em,,_i (s))^ (l - 
which proves the claim (3.46). Thus, we can iterate in ( |3.46| ) with respect to k > 1 and obtain that 

^ / \ 2^ 7 

(3.49) 


^sup^F;^,(t)< (c( 2 ")') ••■••(C( 2 )' 


(sup||u(s)||2,n) 

S>T 


< C 


:2^Ef=7Aj2(52'=ELi 


(sup||u(s)||2_f2) 

s>r 


Therefore, we can take the 1 + = 2^-th root on both sides of (3.49) and let k —>■ + 00 . Using the 

facts that ( := ^ < 00 , C •= ^ < 00 , we deduce 


sup ||m (t) 11^ ^ < lim sup {Emk (t))^/(^+”"'=) < C^2'^^ sup ||n (s) 

fc-7 + OOt>ig s>T- 


i2,n' 


t>to = T‘ 

This clearly proves the proposition 


□ 
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Combining estimate ( 3.35|) (which is also satisfied by the ”generalized” solution) with Proposition 


3.10, and arguing as in the proof of Theorem 3.4, we obtain the following. 


Corollary 3.11. Let the assumptions of Proposition 3.8 he satisfied. For every uq G Lfi (fl), the 


(unique) orbit u{t) = S (t) uq is a strong solution for t > p, for all p > 0. 


Next, since we already know that there exists an absorbing set 13 C (fl) (cf. Ineq. (3.35)) 


for the dynamical system (5 (t) , (f^)), it will be important to show that the semigroup S is also 

asymptotically smooth. This is essential in the attractor theory and the related results we shall 
present in the next section. Clearly, we want to make use of the fact that the weak solutions are 
sufficiently smooth on the interval [/?, oo), for any p > 0. It will suffice to derive a uniform bound 
in Vq (n,(9fl,//). Since the argument relies on the use of key test functions (more precisely, we 
need to take ^ = dtu{t) into the variational equation (3.2)), we will actually need to require more 
regularity of the strong solution, i.e., u G lTj^’®((0, oo); Vq (11, 914, /r)), s > 1. However, lacking any 
further knowledge on these solutions (note that H is an arbitrary bounded open set), we will work 
with ’’truncated” solutions which can be obtained with the Galerkin approach and has independent 
interest. 


3.8 


Lemma 3.12. Let the assumptions of Proposition 
u{t) = s (t) Uq of satisfies 

u G L°°{[p, oo); V 0 ^(H, dLl, p)) n ([p, oo)]Lfi (H)) 
for every p > 0, and the following estimate holds: 

sup + ^ <6*5, 

for some constant Cs = C (p) > 0, independent of t and initial data. 


be satisfied. Then, for uq G (H) any orbit 


(3.50) 


Proof. We recall that by Theorem 


2.15 


is a positive and self-adjoint operator in Lp‘{Ll). 


Then, we have a complete system of eigenfunctions for the operator in Lfi (H) with 

G D{Aq,^^) C V 0 ^(H, p), and G L°° (H). According to the general spectral theory, the 

eigenvalues A* can be increasingly ordered and counted according to their multiplicities in order 
to form a real divergent sequence. Moreover, the respective eigenvectors turn out to form an 
orthogonal basis in V 0 {Ll,dQ.,p) and Lfi (H), respectively. The eigenvectors may be assumed 
to be normalized in Lfi (H). To this end, we can now define the finite-dimensional spaces 

Vn = span {^1, Cn} > 'Poo = 'Jn=lPn- 

Clearly, Voo is a dense subspace of V 0 ^(H, p). As usual, for any n G N, we look for functions 

of the form 

n 

Unit) = ^ei{t)ii (3.51) 

i=l 

solving a suitable approximating problem. More precisely, for any given n > 1 we look for C^-real 
functions e* (•), i = 1,..., n, only depending on time, which solve the approximating problem P (re) 
given by 

(dtUn, 'f')2,0 “b (^©.At'^n; ' 1 ^) 2 ,D + if i^n) 1 'I')2,D ~ *-), (3.52) 

with initial condition 

(re„ (0), 41)2 Q = (re„0)'l')2 Q , (3.53) 
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for all 'if G Vn- Here Uno is the orthogonal projection of uq onto Vn- Observe that 

lim Uno = uq, in (fi). (3.54) 

n^oo 

Notice that by the Canchy-Lipschitz theorem, one can find a unique maximal solution 

UnGC\[0,Tny,D{Ae,f,)nL°° (H)) 


to ( 

3.52 

)-( 

3.53 

), 

for some Tn G (0,r). As in the case when is smooth (see the monographs 



), 

the existence of a generalized solution, defined on the whole interval [ 0 , T], for 


every n G N, can be also obtained with this approach. In particular, notice that (3.35) is also 

satisfied by Un- Here we are interested to derive the bound (3.50). Notice that the key choice of 

1 , 2 / 


the test function ^ = dtUn G ^([O, T]; V 0 (0, cffl, ^u) n (H)) into the variational equation (3.2 ) is 
now allowed for these truncated solutions. We deduce 

^ (t) llvQ’^(r2,9n,/.i) ’ ^^ 2 , 0 ) ^ \\dtUn (f-) Il2,0 ~ 

for all t > 0. Here and below, F denotes the primitive of /, i.e., F (s) = f (y) dy. Multiply this 
equation by t > p > 0 and integrate over (0, t) to get 

t (ll«n Wllv4’2(n,9f2,/.) + ’ 1 ) 2 , 0 ) + ^ (■s)lllo ds 

= [ (lliin (s)llvL 2 (f^ + {F {un (s)), 1 ) 2 ,n) ds, 

for all t > p. Recalling that, due to (H2)-(H3), F is bounded from below, independently of n, and 
\F (s)| < C (1 + |s|^), we infer from ( |3.35 ), 

+ 2 [ < c , (3.55) 


\Un (f^) Ik 


iv4’Ro,ao,/.) 

for some constant c > 0 independent of t, n and p. On the basis of a lower-semicontinuity argument, 
we easily obtain the desired estimate (3.50). The proof is finished. □ 


Notice that the uniqueness of the weak solutions (see Proposition 3.8) is important here. We 
will eventually need to use the fact that each such solution u belongs to L°°{[p,oo)-, L°° (Q) n 
V 0 ^(n, cIH, p)), for every p > 0 , in order to apply a sequence of abstracts results in the next section 
and obtain the desired ’’finite-dimensional” description of the long-term dynamics of (!.!)-(1.3) 
in terms of attractors which possess good properties. Indeed, even though the weak solution of 


Dehnition 3.1 can be also constructed with aid from the Galerkin approach presented above, the 


application of this scheme seems quite problematic for the proof of Proposition 3.10| (see also the 
proof of Theorem 3.4) and it cannot be longer used in that context. In other words, the procedure 


of approximating by the strong solutions Ug gives a weak solution u, while the usual limit procedure 
in the Galerkin truncations Un gives another weak solution u. The solution u coincides with u, 


if uniqueness is known for (1.1)-(1.3). Here is where the two theories seem to depart if the latter 
property is not known. 


4. Finite dimensional attractors 


The present section is focused on the long-term analysis of problem (1.1)-(1.3). We proceed to 
investigate the asymptotic properties of (1.1)-(1.3), using the notion of a global attractor. We begin 
with the following. 

Definition 4.1. A set GB,fi F L? (H) is a global attractor of the semigroup S (f) on L‘^{Q) associated 
with (1.1)-(1.3) if 
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• is compact in L? (fi); 

• Gq,^ is strictly invariant, that is, S{t)Ge,ij, = Ge,^ V i > 0; 

• Ge,f_i attracts the images of all bounded subsets of (0), namely, for every bounded subset B 
of Lf {Gl) and every neighborhood O ofGe,fi in the topology of (O), there exists a constant 
T = T{B; O) > 0 such that S{t)B C O, for every t >T. 

The next proposition gives the existence of such an attractor. 


Proposition 4.2. Let the assumptions of Proposition 3.8 be satisfied. The semigroup S (t) on 


LP‘{Q.) associated with the reaction-diffusion system (1.1)-(1.3) possesses a global attractor Ge,^ in 


the sense of Definition 4-1. As usual, this attractor is generated by all complete bounded trajectories 


of (1.1)-(1.3), that is, Ge,fi = ^|t=oj where JC is the set of all strong solutions u which are defined 
for all t G M-i- and bounded in the (hi) n {Gl, dGl, p)-norm. 

Proof. Due to the dissipative estimates (3.35), (3.45) and (|3.50), the ball 


Bo = £ X := (D) n Ve’^(D,9D,^) : ||n||_:^, < i?} , 


for a sufficiently large radius i? > 0, is an absorbing set for S (t) in (D). Indeed, in light of (3.35), 
it is not difficult to see that, for any bounded set B C (D), there exists a time = t^{B) > 0 
such that S (t) B G Lf (D), for all t > Nex t, we can choose t = t -\- 2 p with t = t,, and /r = 1, 
so that the L?‘-L°° smoothing property (3.45) together with the estimate (3.50) entails the desired 
assertion. Obviously, the ball Bo is compact in the topology of Lf (D). Thus, S (t) possesses a 
compact absorbing set. On the other hand, due to Proposition 3.8 (see also (4.6) below), for every 
fixed t > 0, the map S (t) is continuous on Bo in the L^-topology and, consequently, the existence of 
the global attractor follows now from the classical attractor’s existence theorem (see, e.g. |16jl. □ 


Let us now construct a Lyapunov functional for (!.!)-(1.3). 


Lemma 4.3. Make the assumptions of Corollary 3.11 (or Theorem 3.4). Then the functional 

Ce {u{t)) := ^ h(0llv4'2(n,aD,^) + W)’ 1 ) 2,0 > 

has along the strong solutions of (H the derivative 


d f 

-—CQiu{t)) = — / \dtuit)\^ dx, a.e. f > 0. 
dt Jn 

In other words, the functional £© is decreasing, and becomes stationary exactly on equilibria u^, 
which are solutions of the system: 

— Au f (u) = 0 in Gl, d^uda + (u + 0^ {u))dpL = 0 on dG. (4-1) 

Proof. The proof is a simple calculation which relies essentially on the fact that strong solutions 
are smooth enough, see Definition 3.3 and Remark 3.5 □ 


Corollary 3.11| tog ether with Lemma |4.3| can now be used to study the asymptotic behavior of 
the solutions of (!.!)-(1.3) by means of the LaSalle’s invariance principle (see, e.g., [SUES]). To 
this end, to any (weak) trajectory of Eqns. ([TT])-([g we associate the respective (positive) w-limit 


set 


Ai2 := { 2 / G (D) : 3tn — )> 00 , pn G if (D) such that S (tn) ?/n —^ 2/ in L^-topology} . 
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The following lemma states some basic properties of the w-limit sets (independent of any Lya¬ 
punov function). Its proof is immediate owing to the continuity properties of the strong solu¬ 


tions (see Definition 3.3, Corollary 3.11 and Remark 3.5) and the compactness of the embedding 

V0’^(D,aD,//) -A (j]). 

Lemma 4.4. (i) Any oj-limit set is nonempty, compact and connected. 

(ii) The trajectory approaches its own limit set in the norm of (D), i.e., 

^lim distL2(^Q) {S (t) uq, A+a) = 0. 

(iii) Any co-limit set is invariant: new trajectories which start at some point in A^a remain in 
^-^2 for all times t > 0. 

The dynamical system (5 (t), (D)) is a ’’gradient” system, namely, we have the following (see, 

e.g., [HU Theorem 10.13]). 


Theorem 4.5. Let the assumptions of Proposition 3.8 he satisfied. The global attractor 


coincides with the unstable set of equilibria E^, which consists of solutions of (f.l). 


Thus, the asymptotic behavior of solutions of (1.1)-(1.3) is properly described by the existence 
of the global attractor However, the global attractor does not provide an actual control of 

the convergence rate of trajectories and might be unstable with respect to perturbations. A more 
suitable object to have an effective control on the longtime dynamics is the exponential attractor 
(see e.g. m ). Contrary to the global one, the exponential attractor is not unique (thus, in some 
sense, is an artificial object), and it is only semi-invariant. However, it has the advantage of being 
stable with respect to perturbations, and it provides an exponential convergence rate which can 
be explicitly computed. Our construction of an exponential attractor is based on the following 
abstract result [Ml Proposition 4.1]. 


Proposition 4.6. Let TLjVjVi be Banach spaces such that the embedding Vi V is compact. Let 
M be a closed bounded subset of TL and /et § : B —>■ B be a map. Assume also that there exists a 
uniformly Lipschitz continuous map T : B —>■ Vi, i.e., 

\\Tbi-Tb2\\v,<L\\b^-b2\\n, V6i,52GB, (4.2) 

for some L >0, such that 

||§ 6 i-S 62 ||^ <7 11^1-^211^+ R||T 6 i-T 62 ||v, V 6 i, 62 GB, (4.3) 

for some constant 0 < 7 < ^ and K > 0. Then, there exists a (discrete) exponential attractor 
Aid C B 0 / the semigroup {S(n) := S"^,n G Z_|_} with discrete time in the phase space TL, which 
satisfies the following properties: 

• semi-invariance: § {Aid) C Aid; 

• compactness: Aid is compact in TL; 

• exponential attraction: dist-^(S’^B, Al^) < for aii n G N and for some a > 0 and 

C >0, where distn denotes the standard Hausdorff semidistance between sets in TL; 

• finite-dimensionality: Aid has finite fractal dimension in TL. 


Moreover, the constants C and a, and the fractal dimension of Aid can be explicitly expressed 
in terms of L, K, 7 , ||B||.^ and Kolmogorov’s K-entropy of the compact embedding Vi V, for 
some K = K (L, AT, 7 ). We recall that the Kolmogorov K-entropy of the compact embedding Vi V 
is the logarithm of the minimum number of balls of radius k in V necessary to cover the unit ball 
of Vi (see, e.g., [E]). 

We are now ready to state and prove 










REACTION-DIFFUSION EQUATIONS 


31 


Theorem 4.7. Assume that the nonlinearity f obeys (H2)-(H3) and [i satisfies (H^). The semi¬ 
group S {t) on associated with possesses an exponential attractor in the 

following sense: 

• is hounded in L°° (Q) n V 0 ^(O, dQ, p) and compact in L? (fi); 

• £e,fi is semi-invariant: S (t) ^ 

• £s,ij, attracts the images of bounded (in Lfi {fiX)) subsets exponentially in the metric of L°° (Q), 
i.e. there exist (3 > 0 and a monotone function Q such that, for every hounded set B C 

(Q), 

distL^^n) {S {t)B,£e^^) < Q (||B||i 2 (Q)) e~^\ for all t > 0; 

• £e,^ has the finite fractal dimension in L°° (£1). 

Since an exponential attractor always contains the global one, the theorem implies, in particular, 
the following result. 


Corollary 4.8. The fractal dimension of the global attractor Gq,ii of Proposition 4-^ is finite in 
L°° (n). Moreover, we have 

lim distL^ffi) {S (t) B, Gq ^) = 0, 

i^+oo 

for every bounded set B <Z Lfi (£1). 

We first recall that, due to the proof of Proposition |4.2[ the semigroup S (t) possesses an absorbing 
ball Bo in the phase space X. Thus, it suffices to construct the exponential attractor for the 
restriction of this semigroup on Bq only. In order to apply Proposition |4.6| to our situation, we need 
to verify the proper estimate for the difference of solutions, which is done in the following lemma. 


Lemma 4.9. Let the assumptions of Theorem 
(1.1)-(1.3) such that Ui (0) G Bq- Then the following estimates are valid: 

12 ^—ujt 


Ui (t) - U2 {t)\\l f^ <M\\ui (0) - 1 X 2 (0)|||f^e ‘^^ + K\\ui - , 


hold, and let ui and U 2 be two weak solutions of 

(4.4) 


and 


\\dtui - dtU2\\^2(^^oMv^^fin,dn,f.)r) + / («) " ^2 ^ (0) - 1 x 2 (0)||2,o > 

^ 0 

(4.5) 

for some cj, zx > 0, M, K,C >0, all independent oft and Ui. 

Proof. The injection V^(n, cin,/x) (If) is compact and continuous. The application of Gron- 

wall’s inequality in (3.43) entails the desired estimate (4.4). The second term on the left-hand side 
of (4.5) can be easily controlled by integration over (0,t) in (3.43). More precisely, we have 


\ui (t) - U2 (t)||2,o + / ||ui (s) - U2 (s)llv4.2(o,an,M) - 11^1 “ “2 (0)|||n 


Furthermore, in light of Corollary 3.11 and estimates (3.45), (3.12), recall that we have 


sup ||XXj 
t>o 


X )0 


L°°{n)r\V^^{n,dQ,ii) 


(4.6) 


(4.7) 


1 2 

for some positive constant C independent of t and xxj. Thus, for any test function ^ G V 0 {£l, d£l, /x), 
using the variational identity (3.2) (which actually holds pointwise for t > 0), we have for the 
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function u ■.= ui — U 2 -, 

{dtu (t),0 = {u (t) ,0- if (ui {t)) - f {u2 (t)),0 

< C* Ik Ikllve’^(n,ar2,/i) ’ 

since / G (M), owing to (4.7). This estimate together with (4.6) gives the desired control on 
the time derivative in (4.5). □ 

The last ingredient we need is the uniform Holder continuity of the time map 11 —)■ 5 (t) uq in the 
L°°-norm, namely, 


Lemma 4.10. Let the assumptions of Theorem J^.l be satisfied. Consider u{t) = S{t)uQ with 
uq G Bq. Then, for every t >0, the following estimate holds: 

||u(t) -u(s)||oo,o < for allt,s > 0, (4.8) 

where \ < 1, C > 0 are independent oft,s, u and the initial data. 


Proof. According to (4.7), the following bound holds for u: 


sup ||n 
t>o 


Consequently, by comparison in (3.2) and from (3.50), we have that 


sup 

t>o 




< C 


which entails the inequality 


\u (t) — u (s )||2 Q < C \t — , for all t,s > 0. 


(4.9) 


Inequality (4.8) is a consequence of (4.9) and the L^-L°° smoothing property (3.45). This follows 
from the fact that the nonlinear function / is continuously differentiable. Indeed, due to the 
boundedness of u{t) G L°° (12), a.e. t > 0 and u{t) G V" 0 ’^(r 2 , cAI, /r), a.e. t > 6 > 0 (cf.. Lemma 


3.12), the nonlinearity / becomes subordinated to the linear part of the equation (1.1) (no matter 
how fast it grows). More precisely, obtaining the L?-L°° continuous dependance estimate for the 
difference u{t) — u (s) of any two strong solutions u{t) ,u (s), corresponding to the same initial 


datum, is actually reduced to the same iteration procedure we used in the proof of Theorem |3.4 
The proof is completed. □ 

We can now finish the proof of the main theorem of this section, using the abstract scheme of 
Proposition |4.6[ 


Proof of Theorem 4^ First, we construct the exponential attractor of the discrete map S (T*) 
on Bq (the above constructed absorbing ball in L°° (12) n V 0 ^(r 2 , dLl, /r)), for a sufficiently large T*. 
Indeed, let Bi = [Ut>'r*5 (t) Bo]j^ 2 , where [■]]^2 denotes the closure in the space (12) and then set 
B := 5 (1) Bi. Thus, B is a semi-invariant compact (for the L^-metric) subset of the phase space 
L^(12) and S (T*) : B —)• B, provided that T* is large enough. Then, we apply Proposition |4.6| on 
the set B with Pi = Lfi (12) and S = 5 (T*), with T* > 0 large enough so that Me~‘^'^* = 7 < 5 
(see (4.4)). Besides, letting 

Vi = L^{[0,T*]-,V^/{n,dn,fi)) n W^’\[0,T*] ; (V^’'(12,ai2,//))*), 

V = L^{[0,T*]-,L^ (12)), 

we have that Vi V is compact, with reference to Theorem |2. 15 Secondly, define T : B —)• Vi to 
be the solving operator for (!.!)-(1.3) on the time interval [0,T*] such that Tuq ■= u G Vi, with 
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u (0) = uo € B- Due to Lemma 4.9, (4.5), we have the global Lipschitz continuity (4.2) ofT fromBto 


Vi, and (4.4) gives us the basic estimate (4.3) for the map S = 5 (T*). Therefore, the assumptions 


of Proposition 4.6 are verified and, consequently, the map S = 5 (T*) possesses an exponential 
attractor Aid on B. In order to construct the exponential attractor for the semigroup S{t) 
with continuous time, we note that, due to Proposition |3.8t this se mig roup is Lipschitz contin uous 
with respect to the initial data in the topology of (D), see also (4.6). Moreover, by Lemma 4.10 
the map (t, uq ) i—)> S (t) uq is also uniformly Holder continuous on [0, T*] X B, where B is endowed 
with the metric topology of (H) (actually, even with respect of the metric topology of L°° (H) 
due to the L^-L°° smoothing property of S and the fact that / G Hence, the desired 

exponential attractor for the continuous semigroup S{t) can be obtained by the standard 
formula 

£e,^.= U S{t)Md. (4.10) 

te[o,T*] 

Finally, the finite-dimensionality of follows from the finite dimensionality of Aid and the L^- 
smoothing property of the semigroup S (t). The remaining properties of £b,^i are 


“ n Ve’^) 


{L' 

immediate. Theorem 4.7 is now proved. 


□ 


We conclude with the following result for the dynamical system {T (t) , (H)), associated with 


strong solutions of problem (1.1)-(1.3) (see Definition |3.3[ ). 

Theorem 4.11. Assume that the nonlinearity f obeys (HI) and /i satisfies (H^). The semigroup 
T (t) on L°° (H) associated with strong solutions of possesses an exponential attractor 

ye,ii the following sense: 

• ye,fj. is bounded in L°° (H) n D (H©^^) and compact in Lfi (H); 

• 3^0 ,is semi-invariant: T {t) (Te,/^) C t > 0; 

• 3^0,/i attracts the images of bounded (in L°° (£1)) subsets exponentially in the metric of 
L°° (H), i.e. there exist /3 > 0 and a monotone function Q such that, for every bounded set 
BcL°° (H), 


distiooi^n) (T(t) B, T 0 .;,) < Q (\\B\\ ^oo(s^) 


-ht 


for all t > 0; 


• 3^0 ,has the finite fractal dimension in L°° (D). 

Proof. It is also not difficult to complete the proof of the theorem, using the abstract scheme of 


Indeed, by estimate (3.16), the se migroup T (t) admi ts a bounded absorbing set 

gives a bounded absorbing 


3.10 


Proposition 4.6 

in (H), which combined with the estimate (3.45) of Proposition 
set in L°° (H). Arguing now as in the proof of Lemma 3.12, by using Galerkin truncations of the 
strong solution, we can easily deduce that the same ball Bq is an absorbing set for T (t) in L°° (H). 
Hence, the above scheme applies and we can once again obtain an exponential attractor with 
the stated properties. This completes the proof. □ 

Remark 4.12. (a) All the results in Sections 3-4 also hold if the nonlocal operator in (2.7) is 

replaced by a more general one. 


{Qf,{u),v):= / K{x-y){u{x)-u{y)){v{x)-v{y))d^i^d^x, 

J Jdnxdn 

with a nonnegative symmetric kernel K (x) = fc (|x|) such that K G (T^, dy). This choice 
renders the nonlo cal t erm dissipative everywhere in the estimates, in particula r, se e the 
proof of Theorem 


3.4 


Clearly, the choice k{r) = r 




s G (0,1), gives (2.7). We 


have restricted our attention to this type of kernels only for the sake of convenience and 
simplicity of presentation. 
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(b) All the results of this article also hold with no changes in all the proofs if O C M" is 
an arbitrary open connected set with finite Lebesgue measure. In the definition of all the 
spaces involved, one has only to replace n (7(0) with n (7c(0) where (7c(0) 

denotes the space of continuous functions on 0 with compact support. We have made the 
presentation with 0 bounded only for the convenience of the reader. 

5. Summary 

In this article, we considered a general (scalar) reaction-diffusion equation on iV-dimensional 
bounded domains 0 with non-smooth boundary 90, with N > 2. Our system captures most 
of the specific and variants of the diffusion models considered in the Introduction and analyzed 
in the literature. Our setting also captures a number of additional models that have not been 
specifically identified or analyzed in the literature. We give a unified analysis of the system using 
tools in nonlinear potential analysis and Sobolev function theory on ’’rough” domains, and then 
use them to obtain the sharpest results. In Section we established our notations and gave some 
basic preliminary results for the operators and spaces appearing in our framework. In Section]^ we 
built some well-posedness results for our nonlinear diffusion model, which included existence results 
(Sections 3.1, 3.2), regularity results (Section 3.3), and uniqueness and stability results (Section 
3.2). In Section]^ we showed the existence of a finite-dimensional global attractor and gave some 
further properties, then we also established the existence of an exponential attractor. In addition 
to establishing a number of technical results for our model in this general setting, the framework 
we developed can recover most of the existing existence, regularity, uniqueness, stability, attractor 
existence and dimension results for the well-known reaction-diffusion equation on smooth domains. 

The present unified analysis can be exploited to extend and establish existence, regularity and ex¬ 
istence of finite dimensional attractor results for other important models based on (weakly) damped 
wave equations, systems of reaction-diffusion equations for a vector it = (ui, .■.,Uk) (k > 2 ), par¬ 
abolic problems with degenerate diffusion, and many others. For instance, our framework requires 
only minor modifications to include reaction-diffusion systems for the vector-valued function it; 
the function spaces become product spaces, and the principal dissipation and ’’smoothing” oper¬ 
ators become block operators on these product spaces, typically with block diagonal form. The 
nonlinearities in these models can be treated in a similar way as in Section (see [T6l Chapter II, 
Section 4]). Furthermore, we remark that one can also easily allow f or ti me-dependent external 
forces h (t), h £ (M; (fl)) , acting on the right-hand side of Eqn. (1.1). Indeed, the existence 

results still hold and one can generalize the notion of global attractor and replace it by the notion 
of pullback attractor, for example. One can still study the set of all complete bounded trajectories, 
that is, trajectories which are bounded for all t G M_|_. All the results that we have presented in this 
paper are still true in that case. We will consider such questions in forthcoming contributions. 


6. Appendix 

To make the paper reasonably self-contained, we include supporting material on regularity results 
for the operator and some basic results from monotone operator theory (see, e.g., [52]). 

Let us recall the following result m Corollary 5.3]. 

Theorem 6.1. Let 12 C be an arbitrary open set with finite measure and assume that p 
satisfies (H^). Then the subgradient —Aq^^) generates a strongly continuous (linear) 

semigroup of contractions on In particular, for every uq G L‘^{VL), the orbit 

u{t) = is the unique strong solution of the first order Cauchy problem 

dtu {t, x) + dipQ^ {u (t, x)) = 0, {t, x) G M+ x 12, u\t=o = in 12, 


( 6 . 1 ) 
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such that 

u G (^([O, oo); L^(ri)) n cx)); L^(0)) and u{t, ■) € D{Aq^^), a.e. on (0,oo). 

Moreover, the (linear) semigroup }^>o non-expansive on L°°{Q) in the sense that 

||g-i^e,,.yo||oo,D < lko|loo,D> (6-2) 

for every t>0 and uq G 


We will now state some results for the non-homogeneous Cauchy problem associated with (6.1). 
The first one is taken from |52( Chapter IV, Theorem 4.3]. 

Theorem 6.2. Let ip : H ^ {—oo, -l-oo] be a proper, convex, and lower-semicontinuous functional 
on the Hilbert space H and set A = dp. Let u be the generalized solution of 


u'{t) + A{u) 3 g{t), tG[0,r], 

U\t=0 = Uq, 


(6.3) 


with g G (0, T; H) and uq € D (A). Then p (u) G L^ (0, T), \/tu'{t) G (0, T; H) and u (t) G 
D (A) for a.e. t £ [0, T] . 

The proof of Theorem |3.4| (Section 3) is based on the following regularity result for the evolution 
problem (6.3) with a generic m-accretive operator A. Since we could not find a proof for it in the 
literature, we choose to include here for the convenience of the reader. The theorem is a more 
general version of [521 Chapter IV, Proposition 3.2] with some modified assumptions. 


Theorem 6.3. Let the assumptions of Theorem\6^ be satisfied. Assume that A = dp is strongly 
accretive in H, that is. A — ojl is accretive for some tv > 0 and, in addition, 

g£L^{[5,oo);H)nW^^H[6,ooy,H), 


for every d > 0. Let u be the generalized solution of (6.3) for uq £ D (A). It follows that u{t) £ 
D{A), t > 5 > 0, and the following estimate holds for some v £ D (^4) ; 


^ j ll'^O “ '^Wh + \\9\\l^{[5,oo)-,H) + 




(6.4) 


for each uq £ D (^4) and t > 6 > 0. The constant C > 0 is independent of t, 6, u, v, uq and g. Here 
AP denotes the minimal section of A. 


Proof. We proceed as follows. We first regularize the solution tt by a sequence of approximate 
solutions Ua with a > 0, in which Ua £ ([0, T] ; H) such that Ua (0) = uq £ D (A), and Ua solves 

the following regularized problem 


u'a {t) + Aa (ua (t)) = ga(t), t£ [0, T] . (6.5) 

Here Aa corresponds to the Yosida approximation of A, and g^ £ ([0,T] \H) is a sequence of 

approximate functions such that, as a —)> 0 "^, 

^ 5 in T“ (0, T; H) n (0, T; H). (6.6) 


Note that by [52l Chapter IV, Proposition 1.8] we have Aa = Pa, a > 0, i.e., Aa is Lipschitz 
continuous. Hence by the stan dard Cauchy-Lipschitz theorem, there exists at least one solution 
Ua £ ([0, T] ; H) to problem (6.5). 
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Step 1. If h > 0, then Ua (t + /i) is a solution of (6.5) with pa (t) replaced by pa {t + h), and the 


accretive estimate on Aa (i.e., {A^wi — AaW 2 ,wi — W 2 )u > oj ||t(;i — W 2 \\h , oj > 0) gives 

1 ^ 

2 dt 

< \\ga {t + h)- Qa' 


“ \\Ua {t + h) - Ua + U \\Ua {t + h) - Uo 


H 


1 


\Ua (t + h) - Uo 
1 


H 


< - CJ \\Ua (t + h) - Ua {t)\\H + - WPa (t + h) - Qa 


UJ 


H 


Thus, for 0 < s < t we get, on account of Gronwall’s inequality (see, e.g., [221 Chapter IV, Lemma 

4.1]), 


\ua {t + h) - Ua < IIWa {s + h) - Ua {s)\\]j 


+ U} ^ f e‘^'^ \\ga{T + h) - ga{T)\\]jdT. 

J s 


Dividing both sides of this inequality by h?, then letting /i —?• O'*", we obtain after standard trans¬ 
formations. 


Ua (t) 


< 


H 


< 


Ua {s) 


Ua {s) 


H 


e-^h-s) ^ 


■i: 


,-^h-T) 


9a (t) 


H 


dr 


H 


+ U} 


-1 


9a (t) 


H 


dr, 


for all t > s > 0. 

Step 2. Fix now v £ H and a > 0, and define 


V’(w') = +aiw) - +aiu) - {AaV,W - v)jj , W £ H. 
Note that if: {w) = ip'a {w) — AaV, w £ H, min^og/;- V’ (w) = ip (v) = 0, and 

u'a {t) + Ip' {Ua {t)) = -AaV + ga{t) , t £ [0, T] . 


Since 


(^^p (Ua) , V — Uo^ < Ip (v) — “Ip (Ua) = —Ip {Ua), it folloWS that 


'0 {Ua {t)) < (u^ (t) + AaV - ga{t),V- Uc 


H 


1 d 


= “2^ lb - Ua (t)||// + iAaV,V- Ua (t)) fj - {pa {t),V- Ua ■ 


Hence, by integration over (0, T) we deduce 


(6.7) 


{^p {Ua (t)) + {AaV,Ua (t) - v)+ (pa {t) ,V - Ua (t))iy) dt 


< ^ (lb - UqWh - lb - Ua {T)\\h) ■ 


( 6 . 8 ) 


This is the energ y est imate. In order to estimate the derivative of Ua, we take the scalar product 
in H of equation (6.5) with Ua, multiply the resulting identity by t > 0, to deduce 

t u'aii) ^ (b {Ua (t)) + {AaV, Ua {t) - v)fj + {pa {t) ,V- Ua (t))//)] (6-9) 

= Ip {Ua {t)) + {AaV, Ua {t) -v)jj + {pa {t) ,V- Ua {t)) jj . 
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Integrating over (0,T) and combining with the energy estimate (6.8), on account of Young’s in¬ 
equality, we obtain 

^ (ll'f^-wo||^ + r2p«?^||^ + r^||fi'a(7')ll^) • (6.10) 

The //-norm of (t) satisfies (6.7) with t = T and s = t, so the left-hand side of (6.10) dominates 


T2 


«« (T) 


H 


— UJ 


-1 




5a ( t ) 


H 


drdt. 


( 6 . 11 ) 


We obtain from (6.10) and ( 6.11| ), for each 0 < S < t < T, 


T2 




< - 
H - 2 


\v — Uo\\\j + ll^at’lllf + \\ga (Y)ll^ +10 ^ — 


2 , 2 
9, 




which coupled together with the basic inequality {a — b)^ > ^— b'^, and equation (6.5), yield 


L2{S,T-,H) 


— ||71„ (Ua (t))||^ < A 11^^ - uorn + + 2r2 \\g^ (r)ll^ + cu-V 


Finally, by virtue of (6.6), we deduce the following bound for Aa (ua) in {6,T;H): 

2 


\Aa {ua {T))\\% < ^\\v- uqWh + 2 + 4 \\g (T)ll^ 2a; 


-1 


9 


L^{5,T-,H) 


( 6 . 12 ) 


(6.13) 


for every a > 0 and T > 6 > 0. Arguing now in a standard way (see the proof of [52l Chapter 
IV, Proposition 3.1]), we can pass to the limit as a —)> 0 in (6.13), along a subsequence an —> 0, 
Aa„ {uan (T)) —^ 2/ £ ^{u{T)), so that (6.4) follows from (6.13) by the weak lower-semicontinuity 
of the norm. This completes the proof of the theorem. □ 

Finally, the following corollary follows as a consequence. 


Corollary 6.4. Let the assumptions of Theorem 6.3 be satisfied. For each uq € D (A), there exists 
a unique solution u of (6.3) such that 

ueL°^ ([d, ooy,D (A)) n ([<5, cx)); H ), 

for every (5 > 0. 
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